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I  INTRODUCTION 


A.  STATEMENT  OF  THE  PROBLEM  STUDIED 

As  stated  in  the  original  proposal,  the  focus  of  this  research  is  to 
"...provide  the  data  and  predictive  equations  necessary..."  for  "...accurate 
diagonostic  measurements  of  temperature,  pressure,  and  species  concentration 
by  CARS  spectroscopy...".  These  data  and  predictive  equations  are  to  be  used 
to  describe  the  Raman  vibrational  Q  branch  for  the  molecules  CO  and  H2  as 
functions  of  temperature,  pressure,  and  composition. 

Within  the  context  of  this  proposal,  the  terms  "predictive  equations" 
refer  both  to  validated  forms  for  the  spectral  distribution  functions  and  to 
a  reliable  predictive  basis  for  the  molecular  parameters  which  specify  the 
distribution  function.  Foremost  among  the  relevant  molecular  parameters  are 
the  broadening  and  shifting  coefficients  for  the  individual  transitions  which 
comprise  the  branch.  Although  there  is  some  first  principle  basis  for  the 
predir  .ive  equations,  the  methodology  is  this  area  typically  (but  not 
exclusively  as  we  demonstrate)  is  empirical. 

The  spectral  distribution  functions  need  to  encompass  the  effects  of 
both  velocity  and  phase  changing  collisions.  The  relative  importance  of 
these  two  types  of  collisions  is  very  different  for  the  two  molecules  of 
interest  in  this  work.  In  particular,  the  study  of  the  hydrogen  diatomic 
molecules  yields  important  tests  of  our  fundamental  understanding  of  the 
role  of  velocity  changing  collisions.  The  correlation  of  the  latter  to  phase 
changing  collisions  also  can  be  revealed  by  such  studies.  Prior  to  our 
work,  there  was  relatively  little  known  about  the  temperature  and  collision 
partner  dependence  of  line  formation  for  the  hydrogens. 

For  the  CO  molecule,  the  dominant  line  formation  process  involves  the 
effects  of  rotationally  inelastic  collisions  which  lead  to  a  pressure 
dependent  narrowing  (collapse)  of  the  Q  branch.  In  contrast  to  the  situation 
for  the  hydrogen  diatomics,  there  is  an  established  framework  for 
description  of  these  spectra.  In  this  case,  the  predictive  basis  involves 
the  use  of  rate  laws  for  rotational  energy  transfer  to  specify  a  relaxation 
matrix  formulation  of  the  Q-brancn  spectral  distribution.  However,  the 
efficacy  of  this  basis  must  be  tested  against  accurate  measurements  as 
functions  of  temperature  and  pressure.  Additionally,  first  principle 
calculations  can  be  used  to  provide  a  physically  realistic  assessment  of  the 
rate  law  models.  These  two  areas  were  the  core  of  our  research  on  the  CO 
system . 

B.  BRIEF  OUTLINE  OF  THE  RESULTS  OF  THIS  STUDY 

In  the  following  section  (II) ,  we  present  summaries  of  the  experimental 
results  for  both  the  CO  and  H2  systems.  These  data  involve  extensive 
measurements  of  the  the  rotational  quantum  number  and  temperature  dependence 
of  line  broadening  and  shifting  coefficients  and  detailed  comparisons  of  the 
predictions  of  various  spectral  models  to  spectral  distribution  functions 
determined  over  large  ranges  of  temperature,  pressure,  and  collision  partner. 
In  the  case  of  the  hydrogens,  previously  unobserved  lineshape  anomalies, 
which  are  pressure,  temperature,  and  concentration  dependent,  are  reported. 
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Because  of  the  different  level  of  our  understanding  of  the  physics  of 
line  formation  for  the  two  systems,  the  comparison  of  these  data  to  the 
predictions  of  theory  is  discussed  from  somewhat  different  points  of  view  in 
Section  III. 

For  the  CO  system,  the  discussion  in  Section  III  focuses  on  the  use  of 
energy-gap  rate  law  models  and  the  simultaneous  prediction  of  line  broadening 
coefficients  and  spectral  collapse.  In  this  study  we  are  interested  in  the 
accuracy  of  our  prediction  of  the  details  of  the  CO  Q-branch  spectrum  as  a 
function  of  P  and  T.  It  is  shown  that  the  familiar,  energy-corrected  sudden 
(ECS)  scaling  approach  is  not  consistent  with  the  data  on  CO.  However,  the 
Raman  spectral  data  are  adequately  fit  by  two  different  exponential  gap  law 
models  which  have  raoically  different  assumptions  about  the  role  of  dipolar 
symmetry  forces  in  the  line  broadening  process.  We  are  able  to  resolve  this 
ambiguity  through  a  first  principle,  semi-classical  line  broadening 
calculation  which  we  describe  in  Section  IV. 

For  the  hydrogens.  Section  III  focuses  on  situations  where  spectral 
models  appear  reliable  and  therefore  we  need  to  predict  certain  molecular 
line  shape  parameters.  The  majority  of  the  discussion  involves  attempts  to 
account  for  the  J  and  T  dependence  of  the  line  broadening  coefficients  in 
terms  of  both  vibrational  dephasing  and  rate- law-modelled  inelastic 
collisions . 

In  Section  IV,  we  turn  from  predictive  equations  to  first  principle 
theory.  The  CO  line  broadening  calculation,  mentioned  above,  is  described. 
This  newly  developed  computational  tool  is  shown  to  give  remarkably  accurate 
agreement  with  infrared  and  Raman  line  broadening  data  over  a  large  range  of 
J  and  T.  The  state-of-the-art  in  the  calculation  of  line  broadening  and 
shifting  in  the  hydrogens  also  is  reviewed  in  the  context  of  our  experimental 
results.  With  some  minor  exceptions,  we  find  that  first  principle 
calculations  are  reliable  only  for  the  case  of  D2:He. 

Foreign  gas  broadening  of  CO  is  addressed  in  Section  V.  Here  the  line 
broadening  calculation,  discussed  in  Section  IV,  is  applied  to  the  Q  branch 
of  CO  in  a  N2  host  from  300  to  1500  K. 

To  assess  the  accuracy  of  measurements  of  temperature  and  pressure 
derived  from  CARS  diagnostic  Q-branch  spectra,  an  extensive  sensitivity 
analysis  which  utilizes  the  N2  Q  branch  as  a  model  system  is  presented  in 
Section  VI.  From  this  analysis  we  estimate  the  errors  in  T  and  P  which 
arise  from  uncertainties  ("noise")  inherent  in  the  measurement  of 
multichannel,  single-shot  CARS  spectra  of  N2 .  The  approach  developed  in  this 
section  establishes  a  useful  framework  for  response  to  questions  concerning 
the  accuracy  of  the  information  derived  from  spectroscopic  mea.surements . 
Further,  this  methodology  allows  an  assessment  of  the  sensitivity  of  such 
mea.-:urements  to  inaccuracies  in  the  spectroscopic  data-base  and  can  assist 
the  experimentalist  in  the  selection  of  measurement  schemes  which  are 
optimal]'  insensitive  to  "noise"  and/or  systematic  data  uncertainties. 

1  -le  final  (Conclusions)  section,  we  discuss  the  limits  of  our 
abili' ^  to  predict  the  Q-branch  spectrum  of  the  simple  diatomic  molecule 
systems  ^^--i^  '.ed  in  this  work.  We  also  suggest  areas  for  future  research  in 
support  <)•  ..aRS  diagnostics. 
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II  SUMMARY  OF  EXPERIMENTAL  STUDIES 

The  experimental  results  reported  here  were  obtained  using  two  different 
high-resolution  Raman  gain/loss  spectrometers:  a  cw- laser  based  inverse  Raman 
system  with  nominally  10  MHz  resolution  and  a  quasi-cw  inverse  Raman 
apparatus,  which  uses  a  pulsed  Raman-pump  laser,  with  an  effective  resolution 
of  45  MHz.  These  systems  are  adequately  described  in  the  experimental 
sections  of  the  publications  previously  reported  for  this  project,  cf.  for 
example.  Refs.  [P2] ,  [MIO] ,  [P3] .  and  [P5] .  We  note  for  completeness  that 
these  experimental  systems  provide  very  high  accuracy  measurements  of  the 
imaginary  part  of  the  third  order  nonlinear  susceptibility.  The  resonant 
part  of  the  CARS  susceptibility  is  fully  determined  from  these  measurements. 
Highlights  of  the  experimental  results  are  summarized  for  each  molecular 
system  in  the  following. 


A.  CO: CO 

1)  LINE  BROADENING  COEFFICIENTS 

In  Fig.  II.  1  we  display  the  J  and  T  dependence  of  the  line  broadening 
coefficients  of  self-broadened  CO.  The  solid  lines  (denoted  theory)  through 
the  data  points  are  predictions  of  a  modified  exponential  gap  rate  law  (MEG) . 
This  MEG  law  expresses  the  rate  for  upward  rotational  transitions  from  state 
J  to  state  J '  as 


Wj 


J 


P  a  • 


1  +  A  E(J)/kT6  ]2 

1  +  A  E(J)/KT  J  * 


exp 


r-  ;g(E(j')  -  E(j))' 

I  kT 


(II. 1) 


with  a,  P,  and  6  being  adjustable  parameters,  A  a  constant  (a  measure  of  the 
duration  of  a  collision),  E(J)  the  energy  of  the  rotational  level  J,  P 
pressure,  T  temperature,  and  k  Boltzmann's  constant.  The  downward  J'-»  J 
transitions  are  found  from  the  condition  of  detailed  balance.  The  value  A=2 
is  assigned  in  Ref.  [P5] .  The  parameters  p  (=1.346)  and  S  (=1.322)  have  been 
refined  relative  to  those  presented  in  Ref.  [P5]  in  the  course  of  the  work 
discussed  in  Ref.  [P3].  The  coefficient  a  is  a  function  of  T  of  the  form 


with  a(To)=  0.01265  cm' ^ /atm  for  Tq=  295  K  and  m=  0.212. 


2)  LINE  SHIFTING  COEFFICIENTS 

Our  measurements  of  line  shifting  were  described  in  Ref.  [P3].  The 
procedures  described  therein  '/ere  repeated  in  the  vicinity  of  600  and  1000  K. 
For  the  determination  of  accurate  values  of  line  shifting  coefficients  it  is 
important  to  utilize  relatively  isolated  transitions  and  to  include  first 
order  corrections  for  line  interference.  This  requirement  results  in  the 
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study  of  the  relatively  weak,  high  J  transitions;  however,  we  were  able  to 
establish  reasonably  accurate  values  and  set  limits  to  the  J-dependence  of 
the  shifts.  Since  the  latter  dependence  is  very  small,  we  can  establish  a 
single,  temperature  dependent  shift  coefficient  for  the  entire  Q  branch. 
Ail  the  results  are  displayed  in  Fig.  II. 2. 

These  experimental  findings  are  corroborated  by  semi-classical 
calculations  of  the  shift  coefficients,  see  Section  IV.  The  calculated  Q- 
branch  shifts,  based  only  on  a  parameter  determined  from  the  room  T  shift, 
are  shown  as  the  solid  line  in  the  figure.  The  variation  with  T  is  described 
by 


Aiz/AP  =  -0.0032  (Tq/T)  cm'Vatm 


(II. 3) 


with  To=  295  K  with  estimated  uncertainty  of  ±0.001  cm' ^ /atm.  We  note  that 
the  results  for  CO,  discussed  above,  are  completely  analogous  to  those 
obtained  for  N2  [Rl]. 

3)  SPECTRAL  DISTRIBUTIONS 

One  of  the  most  difficult  aspects  of  the  experimental  research  is  the 
accurate  determination  of  the  shape  of  the  collapsing  Q  branch  as  a  function 
of  P  and  T.  These  measurements  require  high  accuracy  in  both  the  frequency 
and  intensity  axes  over  large  ranges  of  dye  laser  wavelength.  In  addition,  we 
found  that  the  pump  and  probe  beams  must  be  precisely  co-focussed  at  their 
crossing  point  in  order  to  minimize  spurious  signals  (i.e.  apparent  power 
modulations  of  the  probe  laser  caused  by  the  pump  laser) .  We  ascribe  these 
spurious  signals  to  beam  steering  effects,  essentially  "schlieren  effects", 
which  therefore  depend  on  the  real  part  of  the  third  order  susceptibility. 
Having  reduced  these  systematic  errors  to  a  minimum,  we  estimate  the 
uncertainies  in  our  measurements  of  the  spectral  distributions  to  be 
approximately  at  the  level  of  2-3%  of  the  peak  signal,  cf.  for  example  the 
results  in  Appendix  A. 

We  have  studied  the  collapse  of  the  Q  branch  at  room  temperature  very 
extensively,  incorporating  the  results  of  three  independent  series  of 
measurements.  This  work  is  discussed  in  detail  in  the  manuscript  included  as 
Appendix  A.  Representative  results  from  the  lineshape  studies  above  room 
temperature  are  shown  in  Fig.  II. 3.  In  this  figure,  the  data  (symbols)  are 
compared  to  the  predictions  of  the  MEG  law  (solid  lines),  Eqs.  II . 1  and  II. 2, 
with  the  shifts  from  Eq.  II. 3.  Residuals  are  displayed,  offset  for  clarity, 
below  the  spectral  comparison.  We  will  return  to  consider  the  analysis  of 
these  data  in  Section  III. A,  below. 
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Shift  coefficient  (10  ^cm  Vatm) 


C0:C0  Q-Branch  Lineshift  vs  T 


T  (K) 
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Figure  II. 2 


•uncertainties  in  the  last  digit  are  indicated  by  numbers  in  parentheses 
*’see  Ref.  [P4]  . 

Only  for  self-shifting  in  H2  have  we  made  an  extensive  study  of  both  the 
J  and  T  dependence.  The  results  are  presented  in  Fig.  II. 6  (this  study  will 
be  published  [M8]).  Finally,  we  note  that  the  T  dependence  of  the  Ar- 
shifting  of  the  H2  Q(l)  line  is  reported  as  part  of  the  study,  discussed 
below,  on  anomalous  line  shapes,  see  Appendix  B. 
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Rotational  Quantum  No.  J 
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Figure  II. 4 
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Rotational  quantum  number  J 
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Figure  II. 6 


3)  SPECTRAL  DISTRIBUTIONS 

The  hydrogen  diatomics  have  long  been  of  interest  as  systems  in  which 
one  could  accurately  test  our  understanding  of  the  way  velocity  changing 
collisions  lead  to  a  reduction  of  Doppler  broadening  with  increasing  density, 
i.e.  the  well  known  phenomenon  termed  Dicke  narrowing.  The  state  of  this 
understanding  was  briefly  reviewed  in  a  paper  [PI]  written  as  part  of  the 
work  under  this  research  project.  Our  subsequent  experimental  work  is 
summarized  in  the  following. 

The  D2 :D2  and  D2:He  systems,  preliminary  results  having  been  reported  in 
Ref.  [P2] ,  were  extensively  studied  at  room  T  over  a  large  range  of  pressures 
and  for  both  forward  and  backward  scattering.  This  study  (which  will  be 
published  [MIO])  took  advantage  of  the  fact  that  forward  and  backward 
scattering  geometries  simultaneously  exist  in  the  multi-pass  mirror  cell  used 
for  the  cw- inverse  Raman  system.  (The  last  mentioned  fact  explains  the 
shapes  of  the  spectral  distributions  reported  in  [P2].)  The  most  important 
result  of  this  investigation,  for  the  purpose  of  this  project,  is  that  the 
soft-collision  model  quantitatively  describes  all  the  pressure,  scattering 
wavevector,  and  J-dependence  of  these  spectra.  For  compltLeness  we  include  a 
short  description  of  this  soft-collision  model. 

This  model  is  a  solution  to  a  Fokker-Planck  type  kinetic  equation  for 
the  evolution  of  phase  in  the  presence  of  both  velocity  and  phase  changing 
collisions.  The  frequency  dependence  of  the  third  order  susceptibility  we 
use  is  then 


-k^v^ 

X  «  /dtexp(-i[w-W(,  ]  t)exp(-Ft)exp( - {Tt-l+exp(-Tt) )) 

2t2 

(II. 

•-1 

II 
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H 

(II. 6) 

k  =  1  kp  -k.  1 

(II. 7) 

where  F  and  T  are,  respectively,  the  phase  and  velocity  changing  collision 
rates,  p  is  the  density,  k  is  the  magnitude  of  the  momentum  transfer,  and  v 
is  the  most  probable  speed  =(2kgT/m)^  .  In  Eq.  II. 4  w  is  the  pump-probe 
frequency  difference  and  the  resonance  occurs  at  .  This  soft  collision 
lineshape  has  the  proper  Gaussian  lineshape  characteristic  of  Doppler 
broadening  at  very  low  densities,  i.e.  T  -►  0,  and  at  higher  densities 

approaches  the  Dicke  limit.  The  latter  results  at  densities  such  that  the 
quantity  k  =  k^v^/(2T^)  «  1.  In  this  limit  the  lineshape  is  Lorentzian  with 
width  (HWHM) 

Dgk^ 

Fg  =  —  +  IP  (II. 8) 

with 
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(II. 9) 


Do  = 


the  optical  diffusion  constant.  The  magnitude  of  k  determines  deviations 
from  the  Dicke  limit.  Increasing  k  will  give,  at  lowest  order,  corrections 
which  cause  the  lineshape  to  become  slightly  sub-Lorentzian,  i.e.  corrections 
which  will  allow  the  shape  to  approach  the  Gaussian,  low-density  limit.  k 
depends  inversely  on  the  square  of  the  density  and  directly  on  the  square  of 
the  momentum  transfer.  Variable  density,  forward  and  backward  scattering 
allows  a  careful  test  of  the  accuracy  of  the  lineshape  function. 

An  observation,  not  yet  reported  in  the  literature,  relevant  to  the 
applicability  of  the  soft-collision  model  to  the  D2  lineshape  is  shown  in 
Fig.  II.  7.  The  host  gas  in  this  case  is  Ar,  and  the  principal  scattering 
geometry  is  backward.  The  concentration  (mole  fraction)  of  D2  is  50%  in  this 
approximately  4  atm  room  T  sample.  The  experimental  spectra  and  theoretical 
models  displayed  in  Fig.  II. 7a  show  the  forward  and  backward  components  which 
make  up  the  observed  spectrum  and  the  residual  of  the  fit.  The 
characteristic  shape  of  this  residual  can  be  used  to  infer  that  the  data 
would  be  better  fit  by  a  complex  soft-collision  model,  cf.  [R2].  The  complex 
soft-collision  model  (in  which  T,  Eq.  II. 6,  is  a  complex  number)  results  from 
a  solution  to  the  generalized  Fokker-Planck  equation  for  the  time 
development  of  the  optical  coherence.  The  result  shown  in  Fig.  II. 7b 
indicates  the  improvement  in  accuracy  provided  by  this  spectral  distribution 
function.  The  latter  form  of  the  lineshape  function  is  expected  under 
conditions  where  there  are  correlations  between  velocity  and  phase  changing 
collisions.  Further  evidence  for  such  a  correlation  is  presented  below  when 
we  discuss  anomalous  lineshapes  in  the  pressure  broadened  regime  for  the 
Dj :Ar  and  Hj :Ar  systems. 

Studies  of  the  T  dependence  of  the  lineshapes  of  the  hydrogen  diatomics 
have  focussed  mainly  on  the  Hj :Ar  system;  however  we  have  determined  the 
lineshape  as  a  function  of  P  at  1000  K  for  the  self -broadened  Q(l)  transition 
of  H2 .  The  results  are  consistent  with  the  Dicke  model,  which  is  the  form  of 
the  soft-collision  model  in  the  density  region  studied.  The  lineshape  is  a 
symmetric  Lorentzian  whose  half-width  is  shown  in  Fig.  II. 8.  The  effective 
diffusion  constant  in  this  case  is  not  equal  to  the  mass  diffusion  value,  in 
sharp  contrast  to  all  room  T  observations  on  this  system,  see  Ref.  [PI].  It 
is  difficult  to  interpret  this  result  at  our  current  of  understanding  line 
formation  in  the  hydrogens,  see  the  following  and  Appendix  B. 

Our  studies  of  D2 :Ar  and  H2:Ar  have  revealed  previously  unknown 
characteristics  of  the  spectral  lineshapes  and  their  dependence  on  pressure, 
concentration,  and  temperature.  These  are  discussed  in  detail  in  Appendix 
B.  The  essential  physics  which  explains  these  results  has  two  important 
aspects.  The  first  involves  the  fact  that  for  light,  optically  active 
molecules  in  a  bath  of  heavy  perturbers  the  collision  energy  is  strongly 
correlated  with  the  laboratory  speed  of  the  optically  active  molecule.  Since 
shifting  (and  broadening)  coefficients  are  functions  of  the  collision  energy, 
we  encounter  a  situation  where  there  potentially  is  a  significant 
inhomogeneous  component  in  the  observed  lineshape.  This  is  the  source  of  the 
asymmetric  lineshapes  shown  in  Figs.  II. 9  and  II. 10.  Under  conditions  where 
the  collisional  dephasing  rate  is  small  relative  to  the  speed  changing  rate, 
and  when  the  latter  rate  is  of  the  order  of  the  reciprocal  of  the 
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Figure  11.7a 
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Figure  11.7b 
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Figure  II. 8 
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Figure  11.9 
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18  Figure  II . 10 


inhomgeneous  width,  we  have  the  possibility  of  collisional  narrowing  of  the 
speed  dependent  inhomogeneous  width,  in  analogy  to  Dicke  narrowing,  discussed 
above.  This  second  aspect  of  the  process  of  line  formation  is  evidenced  by 
the  concentration  dependence  of  the  lineshapes,  cf.  Figs.  II. 9  and  II. 10,  and 
the  nonlinear  variation  of  linewidth  with  concentration  shown  in  Fig.  II. 11. 
An  analytical  model  which  incorporates  the  two  aspects,  discussed  above,  and 
which  is  consistent  with  a  large  range  of  experimental  observations  on  the 
H2 ;Ar  system,  is  described  in  Appendix  B. 
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Figure  11.11 
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Ill  COMPARISON  TO  PREDICTIVE  EQUATIONS 


The  state  of  our  knowledge  of  the  Q-branch  spectrum  is  much  more  highly 
developed  in  the  case  of  CO  than  for  the  H2  molecule.  In  the  former  case  the 
broadening  mechanisms  and  spectral  distribution  function  are  reasonably  well 
understood  and  the  research  involves  the  development  of  a  reliable  predictive 
basis  for  the  fundamental  parameters  which  specify  the  spectrum.  For  the 
diatomic  hydrogens  there  are  questions  with  respect  to  the  fundamental 
"pressure  broadening"  mechanisms,  e.g.  rotational  inelasticity  versus  elastic 
vibrational  dephasing,  and  with  respect  to  the  form  of  the  spectral 
distribution  function,  for  example,  what  is  the  spectral  shape  in  the 
presence  of  simultaneous  velocity  and  phase  changing  collisions?  The 
research  reported  here  (Section  II.B.3.b  and  Appendix  B)  has  revealed  another 
contribution  to  the  Hj  lineshape  in  the  "pressure  broadened,  impact  regime". 
In  this  section,  we  will  concentrate  on  an  enumeration  of  the  results  of 
tests  of  our  ability  to  predict  the  Q  branch  spectrum  for  those  cases  where 
the  broadening  mechanisms  and  form  of  the  spectral  distribution  function  is 
reasonably  well  verified. 


A.  CO; CO 

1)  ROOM  TEMPERATURE 

Although  our  first  studies  [P3,  P5]  of  the  CO  system  employed  the  MEG 
law,  we  subsequently  expanded  consideration  to  other,  well-established  rate 
law  models  for  rotationally  inelastic  collisions.  The  details  of  this  work 
are  discussed  in  Appendix  A.  A  conclusion  derived  from  this  work  is  that  the 
energy  corrected  sudden  scaling  law  analysis  (ECS)  of  CO  [R3],  although  it 
gives  good  account  of  the  J-dependence  of  the  line  broadening  coefficients, 
does  not  properly  describe  the  collapse  of  the  Q  branch  with  increasing 
pressure  at  room  temperature.  In  contrast  to  this  result,  we  also 
demonstrated  that  the  modified  exponential  energy  gap  law  (MEG)  and  the 
statistical  power-exponential  energy  gap  law  (SPEC)  cannot  be  distinguished 
on  the  basis  of  agreement  with  both  linewidths  and  Q-branch  collapse.  This 
is  surprising  in  that  these  two  laws  fit  the  data  only  if  radically  different 
assumptions  with  respect  to  the  propensity  for  odd  delta-J  transitions  are 
invoked.  In  the  case  of  the  MEG  law,  the  best  fit  to  the  complete  data  set 
requires  that  all  delta-J  values  are  allowed  with  the  consequence  that  more 
than  50%  of  the  total  width  must  be  associated  with  dipolar  (more  generally 
odd  multipolar)  symmetry  forces.  The  SPEG  law  on  the  other  hand  fits  the 
data  set  best  for  the  assumption  that  only  even  delta-J  transitions  are 
allowed.  Although  the  latter  selection  rule  is  not  required  by  the  symmetry 
of  the  CO  molecule,  the  appearance  of  large  rates  for  odd  delta-J 
transitions  does  appear  to  be  at  variance  with  the  observation  that  CO  is 
almost  homonuclear  in  its  electrical  and  transport  properties.  This 
ambiguity  was  resolved  (in  favor  of  the  SPEG  law  with  quadrupolar  symmetry 
forces)  by  a  first  principle  semi-classical  line  broadening  calculation. 
This  calculation  is  summarized  in  Section  IV. 

The  room  temperature  CO  Q  branch  is  now  well  understood.  The  subtle 
differences  between  the  collapsed  spectra  and  the  SPEG  law  based  relaxation 
matrix  are  at  the  limit  of  the  uncertainties  in  oui  measurement  of  the  Q- 
branch  spectrum.  Further  experimental  verification  of  the  predictions  of 
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theory,  stimmarized  above,  most  likely  will  result  only  from  direct 
measurements  of  state-to-state  rates  for  rotational  energy  transfer.  We  now 
turn  to  extension  of  our  study  of  the  CO  Q  branch  to  elevated  temperatures. 

2)  ABOVE  ROOM  TEMPERATURE 

We  presented  a  reasonably  accurate  MEG  law  based  parametrization  of  the 
J  and  T  dependence  of  the  line  broadening  coefficients  in  Section  II.  In 
light  of  the  discussion  above,  it  appears  that  this  is  a  rather  "unphysical" 
description  of  the  Q  branch.  We  persist  in  the  use  of  the  MEG  law  two 
reasons.  First,  the  description  of  the  line  broadening  coefficients  is 
simple  and  accurate.  A  similar  quality  description  based  on  the  SPEG  law 
does  not  appear  possible.  Second,  the  description  of  the  Q-branch  collapse 
at  elevated  temperatures  is  not  significantly  improved  by  use  of  the  SPEG 
law.  Thus,  our  assessment  is  that  the  MEG  law,  albeit  microscopically 
inaccurate,  does  provide  a  remarkably  simple,  accurate,  and  useful  method  of 
predicting  the  CO  Q-branch  spectrvim  over  a  large  range  of  T  and  P. 

Judging  solely  on  the  basis  of  the  room  temperature  results  we  might 
expect  that  the  SPEG  law  would  be  more  accurate  in  the  description  of  the 
high  T  Q  branch.  We  have  tested  this  idea  by  adjusting  the  SPEG  law  to 
agree  with  the  1050  K  line  broadening  coefficients  (calculated  from  the  MEG 
law  of  Section  II)  and  then  comparing  to  the  experimental  data  on  the  Q- 
branch  collapse  at  =1050  K.  This  result  is  shown  in  Fig.  III.  la.  The 
systematic  differences  between  the  SPEG  model  and  the  experimental  data  are 
quite  comparable  to  those  found  for  the  MEG  law,  cf.  Fig.  II. 3b.  We  observe 
that  the  differences  seen  in  these  two  figures  are  characteristic  of  an 
underestimate  of  the  degree  of  collapse.  Again,  as  with  the  room  temperature 
analysis,  the  two  laws  make  very  different  assumptions  about  the  propensity 
for  odd  delta-J  transitions.  In  an  attempt  to  better  understand  these 
observations ,  we  again  turn  to  the  theoretical  calculations . 

From  the  calculations,  discussed  in  Section  IV,  we  know  that  the  dipolar 
symmetry  forces  make  a  finite  contribution  to  the  total  line  width  at  all 
temperatures.  Further,  we  know  from  the  temperature  dependence  of  either  the 
MEG  or  SPEG  laws  that  the  total  inelastic  rate  spreads  over  a  larger  range  of 
delta-J  as  T  Increases.  Thus,  the  relative  contribution  to  the  collapse  of 
the  Q  branch  from  the  delta-J=±2  and  delta-J=±4  (and  higher)  transitions 
decreases  with  increasing  T.  This  suggests  that  the  small  dipolar  symmetry 
component  in  the  inelastic  rate  may  play  a  more  significant  role  in  the 
collapse  of  the  Q  branch  at  elevated  temperatures.  To  test  this  idea,  we 
examine  the  predictions  of  collapse  for  the  case  where  we  allow  approximately 
5%  of  the  total  linewidth  to  be  associated  with  delta-J=±l  transitions  and 
readjust  the  SPEG  law  (even  delta-J)  to  the  experimental  linewidths  (from  the 
MEG  law,  less  5%)  at  1050  K.  Thus,  we  have  a  relaxation  matrix  which  has  the 
even  delta-J  SPEG  law  symmetry  with  the  addition  of  nearest-neighbor  (delta- 
J=±l)  line  interference  (from  dipolar  symmetry  forces)  and  the  same  total 
inelastic  rate  (linewidth)  as  the  MEG  and  the  SPEG  discussed  previously.  The 
comparison  of  this  "modified  SPEG  model"  to  the  experimental  data  is  shown  in 
Fig.  III. lb.  We  suggest  that  the  improved  agreement  displayed  in  this  figure 
reflects  a  more  physically  realistic  treatment  of  the  microscopic  details  of 
the  collapse.  It  should  be  noted  however  that  this  simple  modification  of 
the  relaxation  matrix  is  not  a  highly  accurate  approximation  for  the  J- 
dependence  of  the  inelastic  contributions  of  the  dipolar  symmetry  forces. 
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V 


The  extension  of  this  analysis  to  a  simple .  accurate  predictive  basis  for  all 
temperatures  does  not  appear  to  be  straightforward  or  particularly  valuable 
in  light  of  the  current  needs  for  accurate  data.  The  current  status  of  our 
understanding  of  the  high  temperature  Q-branch  spectrum  of  simple  diatomics 
will  be  the  subject  of  a  future  publication  [M13]. 


B.  DIATOMIC  HYDROGEN  SYSTEMS 

1)  ROOM  TEMPERATURE 

a)  LINE  BROADENING  COEFFICIENTS 

We  have  shown  that  the  MEG  law  can  provide  a  physically  reasonable  basis 
for  the  description  of  the  J-dependence  of  line  broadening  coefficients  for 
the  D2 :X  system,  see  Section  II  and  Ref.  [P2] ,  and  for  the  HD: HD  system  [P4] . 
In  both  systems,  all  available  linewidth  (both  infrared  and  Raman  for 
multiple  branches)  and  inelastic  rate  data  were  included  in  the  analysis. 
Most  significant  among  the  results  of  these  analyses  was  the  assignment  of  a 
reasonable  value  for  the  contribution  of  (elastic)  vibrational  dephasing  and 
the  identification  of  a  significant  vibrational  state  dependence  in  the 
inelastic  rates  for  HD.  A  similar  parametrization,  utilizing  the  MEG  law, 
was  made  for  our  measurements  of  H2:H2;  however,  the  form  of  the  law  had  to 
be  significantly  modified  (with  no  clear  physical  motivation)  in  order  to 
account  for  the  T-dependence .  This  has  led  to  a  re-examination  of  the  use  of 
rate  laws  for  this  system  which  we  will  discuss  in  Section  III. B. 2. a,  below. 


b)  SPECTRAL  DISTRIBUTION 

For  all  the  hydrogen  diatomics  Doppler  broadening  and  Dicke  narrowing 
are  important  components  of  the  Q-branch  spectral  distribution  function.  For 
the  systems  H2  :H2  [P1,R4]  and  D2  :X  (X=D2 ,  He,  H2 )  [P2,  MIO]  at  room 
temperature,  the  soft-collision  (Galatry)  model  (R5,  R2]  has  been  shown  to 
quantitatively  describe  the  spectrum  with  the  coilisional  narrowing  parameter 
determined  from  the  mass  diffusion  constant.  In  the  case  of  HD: HD,  the  large 
amount  of  pressure  broadening  to  some  extent  masks  the  details  of  the  Dicke 
narrowing  process;  however,  we  were  able  to  quantitatively  predict  the 
spectrum  with  an  optical  diffusion  constant  only  slightly  smaller  than  the 
mass  diffusion  coefficient  [P4] .  To  interpre  this  observation,  it  is 
important  to  remember  that  it  is  the  collision  controlled  diffusion  of  the 
optical  coherence  (created  by  the  applied  fields)  which  we  measure 
spectroscopically,  not  strictly  the  diffusion  in  space  of  the  molecule. 

Our  studies  of  D2  in  collision  with  Ar  (see  Section  II) ,  N2 ,  and  CH^  and 
the  results  for  H2:Ar,  discussed  in  Section  II  and  Appendix  B,  reveal  that  we 
currently  have  no  general  predictive  basis  (for  the  lineshape  or  for  some  of 
the  relevant  spectral  parameters)  for  these  systems. 


2)  ABOVE  ROOM  TEMPERATURE 
a)  LINE  BROADENING  COEFFICIENTS 
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As  suggested  above,  it  was  the  J  and  T-deoendence  of  the  line  broadening 
coefficients  which  motivated  a  re-examination  of  rate  law  models  for  the 
H2  '.^2  system.  Additional  interactions  which  should  be  considered  in  a  model 
for  self -broadening  of  Hj  are  the  strong  role  of  resonant  R-R  transfer  [R6] 
and  the  possibility  of  resonant  (J-dependent)  vibrational  dephasing  [R"^] 
(which  plays  a  strong  role  in  the  shifting  coefficients,  see  Section  IV). 

To  explicitly  allow  for  both  R-T  and  R-R  rotationally  inelastic 
processes,  we  have  invoked  the  use  of  the  ECS  scaling  law  (discussed  and 
referenced  in  Appendix  A) .  Additionally,  a  J-dependent  vibrational  dephasing 
component  was  considered.  The  outcome  of  the  analysis  [M9]  is  shown  in  Figs. 
II. 5  and  III. 2.  The  results  displayed  in  these  figures  are  for  the  case  of 
J-independent  vibrational  dephasing.  The  last  mentioned  model  is  more 
consistent  with  other  (pure  rotational  Raman)  linewidth  data.  This  J- 
independent  vibrational  dephasing  for  the  room  temperature  and  above  data  can 
be  reconciled  with  the  observation  of  resonant  vibrational  dephasing  at  low 
temperature  [R7]  and  is  in  accord  with  the  previously  mentioned  analyses  of 
the  other  hydrogen  diatomics. 

Although,  this  ECS  based  scaling  law  model  is  predictive,  it  is  far  from 
simple.  It  Involves  the  assumption  of  an  exponential  gap  base-rate  for  the 
R-T  process  and  a  single  base-rate  for  R-R  processes.  Each  process  has  a 
characteristic  scaling  length  and  magnitude,  and  for  the  R-T  process  there  is 
an  energy  gap  parameter.  It  is  quite  satisfying,  to  our  physical  intuition, 
that  the  scaling  lengths  and  the  energy  gap  parameters  are  found  to  be 
temperature  independent  and  that  the  two  magnitudes  have  reasonably  simple  T- 
dependent  scaling.  These  results  appear  to  give  a  physically  realiscic  and 
reliable  description  of  the  J  and  T-dependence  of  the  line  broadening  in 
H2 :H2  . 


We  observe  that  modelling  the  various  systems  described  above  required  a 
great  deal  of  information  on  the  specific  system  (state-to-state  rates, 
multiple  spectroscopic  branches).  Wlien  temperature  and  collision  partner 
dependence  are  added  to  the  problem,  we  are  forced  to  conclude  that  our  level 
of  understanding  and  ability  to  predict  are  rather  limited  at  present. 

b)  SPECTRAL  DISTRIBUTION 

The  study  of  Dicke  narrowing  in  self-broadened  H2  at  1000  K,  discussed 
in  Section  II,  illustrates  our  inability  to  quantitatively  predict  even  this 
component  of  the  spectrum  of  the  self-broadened  hydrogen  diatomics.  In  this 
case,  the  optical  diffusion  constant  is  found  to  be  larger  than  the  mass 
diffusion  constant.  This  observation  is  exactly  the  inverse  of  the 
conclusions  of  studies  of  all  other  systems  for  which  a  reasonably  complete 
analysis  has  been  completed.  The  latter  systems  give  results  consistent 
with  our  expectation  that  the  optical  coherence  moves  through  the  gas  as 
fast  as  an  individual  molecule.  The  H2  result  suggests  that  resonant  v-v 
(vibration  to  vibration)  transfer  may  play  a  significant  role  in  enhancing 
the  optical  diffusion  constant. 

The  lineshape  functions  associated  with  the  speed  dependence  of 
broadening  and  shifting  extend  the  need  to  predict  additional  "kinetic" 
coefficients,  e.g.  speed  changing  rates,  even  when  we  are  interested  in  the 
normal  pressure  broadened  regime.  Currently,  there  is  no  predictive  basis 
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Pure  Dephasing 


for  the  complete  range  of  kinetic  and  line  broadening  parameters  necessary  to 
specify  the  spectra  of  the  hydrogen  diatomics  in  an  arbitrary  host  gas  over 
an  extended  range  of  temperature  and  pressure . 
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IV  COMPARISONS  TO  FIRST  PRINCIPLE  THEORY 


The  predictive  equations  discussed  in  the  preceeding  section  are 
essentially  empirical  expressions  having  physically  realistic  bases  but  no 
simple  expression  derived  from  or  value  determined  from  first  principle 
theory.  In  this  section  we  turn  to  a  consideration  of  the  level  of  our 
ability  to  understand  and  predict  the  lineshape  and  its  relevant  parameters 
from  fundamental  theory.  Our  discussion  is  limited  to  a  few  topics 
specifically  relate  1  to  the  molecular  systems  of  interest  in  this  work. 

A.  CO; CO 

1)  LINE  BROADENING 

It  became  clear  early  in  our  research  that  we  would  be  well  served  to 
initiate  a  first  principle  computational  study  of  the  CO  system  to  support 
the  measurement  program.  The  motivation  involved  two  primary  elements. 
First,  there  were  insufficient  experimental  data,  for  example  on  state-to- 
state  rates,  to  allow  a  critical  analysis  of  rate  law  models.  This  situation 
would  become  more  difficult  as  we  extended  measurements  above  room 
temperature.  Second,  the  experimental  difficulties  associated  with  making 
accurate  spectral  measurements  under  conditions  of  foreign  gas  broadening 
would  probably  preclude  our  making  substantial  new  contributions  in  this 
area.  Reliable  theory  could  then  significantly  improve  our  understanding  and 
supply  needed  data. 

Our  approach  to  the  calculation  of  molecular  line  broadening  and 
shifting  was  guided  by  a  few  simple  criteria.  First,  the  theorertical 
formulation  of  the  problem  must  result  in  a  manageable  computational 
requirement  when  applied  to  moderately  heavy  diatom-diatom  systems.  Second, 
this  formulation  must  be  based  on  intermolecular  potentials  which  have  been 
independently  verified  as  accurate  and  which  can  be  readily  incorporated 
into  the  calculational  scheme.  Finally,  the  approach  must  provide  line 
broadening  values  such  that  we  can  test  against  the  wide  variety  of 
experimental  data  available  for  CO  over  a  temperature  range  from  100  K  to 
1500  K.  We  choose  as  our  theoretical  basis  the  well  established  Robert  and 
Bonamy  [R8]  semi-classical  formulation  of  molecular  line  broadening.  The 
salient  features  of  this  formulation  are: 

(1)  through  the  use  of  cumulant  or  linked-cluster  expansion  techniques, 
all  terms  which  are  bilinear  in  the  anisotropic  part  of  the 
intermolecular  potential  are  accounted  for  to  infinite  order  of 
perturbation  theory.  This  yields  a  form  for  the  broadening  cross 

.  section  which  remains  analytic  for  all  impact  parameters  and  therefore 
avoids  the  need  for  any  cut-off  procedures. 

(2)  a  parabolic  trajectory  model,  based  upon  the  isotropic  part  of  the 
intermolecular  potential,  is  used  to  represent  collision  trajectories 
near  the  distance  of  closest  approach. 

(3)  Iciig  range  electrostatic  forces  are  accounted  for  through  the  use  of 
point  multipoles;  in  the  case  of  CO  this  involves  the  dipole 
(essentially  negligible)  and  quadrupole  moments. 
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(4)  the  use  of  a  pairwise,  site-site  Lennard- Jones  6-12  potential  to 
realistically  model  interactions  at  short  range. 


We  note  that  all  terms  given  in  Refs.  [R8-R11]  are  explicitly  included  in 
these  calculations  and  that  the  forms  for  the  resonance  functions  were 
adapted  from  [Rll],  Relavent  molecular  interaction  parameters  are  given  in 
Table  II  of  Appendix  A.  Full  details  of  these  calculations  will  be  presented 
elsewhere  [Mil ] . 

The  results  of  our  calculations  are  compared  in  Fig.  IV.  la  with  the 
experimental  results  of  Nakazawa  and  Tanaka  [R12,  R13]  for  the  infrared  R- 
branch.  In  Fig.  IV. lb,  the  calculations  are  compared  with  our  experimental 
data  for  the  Raman  Q  branch.  From  these  figures,  we  see  that  there  is 
excellent  agreement  between  observed  and  calculated  linewidths  over  a  large 
range  of  J  and  T  for  both  infrared  and  Raman  transitions.  The  temperature 
range  of  these  calculations  involves  over  an  order  of  magnitude  in  mean 
collision  energy,  thereby  probing  a  large  extent  of  the  CO: CO  intermolecular 
potential.  It  is  important  to  note  that  this  potential  is  derived  from  the 
literature  and  not  adjusted  to  predict  the  experimental  linewidths. 

The  quality  of  the  agreement  shown  in  Fig.  IV.  1  is  indicative  of  the 
accuracy  of  the  calculational  tool  developed  for  this  work.  Some  of  the 
critical  features  of  this  tool  are  the  inclusion  of  all  terms  in  the  fourth- 
order  expansion  of  the  atom-atom  potential  and  the  use  of  a  very  fine  grid  to 
evaluate  the  integrals  over  impact  parameter.  We  note  that  the  short-range 
portions  of  the  potential  make  significant  contributions  to  the  total  line 
width  in  CO: CO  and  thus  this  computational  approach  must  provide  a  reasonably 
accurate  representation  of  the  short-range,  anisotropic  potential  and  of  the 
path  for  close  collisions.  These  observations  are  analyzed  and  discussed  in 
more  detail  in  the  forthcoming  publication  [Mil]. 

2)  LINE  SHIFTING 

The  isotropic  Raman  Q  branch  lines,  in  contrast  to  other  spectroscopic 
branches,  do  not  exhibit  elastic  reorientational  phase  shifts.  Line  shifts 
for  these  transitions  arise  from  the  vibrational  state  dependence  of  the 
isotropic  part  of  the  intermolecular  potential,  and  in  our  analysis  we  treat 
this  in  the  lowest  order  (linear  in  ^  act  ion^  perturbation  theory. 
As  described  in  Appendix  A,  the  isotropic  potenial  is  represented  by  a 
Lennard-Jones  model.  Line  shifting  is  accounted  for  via  a  vibrational  state 
dependence  of  the  Lennard-Jones  well-depth  e.  A  fractional  change, 
(Se/f )=0 . 55% ,  is  found  to  give  agreement  with  the  experimentally  determined 
shifts  at  295  K.  From  the  calculation  we  verify  the  experimental  observation 
that  the  shift  is  very  weakly  J-dependent  (the  calculated  values  vary  less 
than  0.001  cra‘^/atm  variation  from  J=0  to  J=24) .  Further,  this  vibrational 
dephasing  makes  less  than  0.0005  cm" ^ /atm  contributions  to  the  linewidths. 
Finally,  the  calculation  predicts  a  T' ^  scaling  of  the  shift  with 
temperature.  These  observations  are  in  agreement  with  all  the  experimental 
results  presented  in  Section  II. A. 2. 


B.  DIATOMIC  HYDROGEN  SYSTEMS 
1 )  LINE  BROADENING 
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Although  the  H2  molecule  was  the  subject  of  the  very  first  calculation 
of  a  Raman  Q-branch  linewidth  [R6] ,  there  remain  many  gaps  in  our  ability  to 
predict  line  broadening  in  the  hydrogen  diatomic  systems.  We  have  reviewed 
much  of  the  theoretical  work  on  these  systems  and  made  comparisons  to  our 
experimental  studies  in  Refs.  {P2]  and  [P^]  .  The  most  complete  and 
successful  first  principle  calculation  has  recently  been  made  to  account  for 
our  measurements  on  D2:He.  The  theoretical  work  grew  out  of  discussions  with 
Dr.  Louis  Monchick  of  the  Applied  Physics  Laboratory  of  The  Johns  Hopkins 
University.  This  work  has  resulted  in  three  separate  publications  which 
include  line  broadening  [R14,  R15] ,  line  shifting  [R15],  and  the  shape  of  the 
Dicke-narrowed  spectrum  [R16] . 

The  calculations  of  Monchick  and  co-workers  are  based  on  a  completely  ab 
initio  potential  surface,  use  "standard"  (S-matrix  based)  expressions  for  the 
cross-sections,  and  determine  the  latter  from  converged,  quantum  mechanical, 
close-coupling  scattering  calculations.  The  computational  dimensions  of  this 
problem  are  tractable  because  we  are  interested  in  D2  ,  with  only  a  few 
populated  rotational  levels  at  room  temperature,  in  collision  with  a 
monatomic  equal  mass  host. 

In  their  second  publication  [R15]  they  provide  the  following  table  (we 
add  our  vibrational  dephasing  estimate  and  slightly  change  the  experimental 
shifts  for  J=3  and  4,  cf.  Table  II. 1): 


Table  IV. 1 

Line  broadening  and  shifting  coefficients  for  D2 :He  at  295  K 
(units  10'^  cm'^/amagat) 


7  exp't  7  theory  shift 


inelastic 


J 

total 

deph. 

total 

v=0 

V=1 

deph . 

exp ' t . 

thy. 

0 

2.35 

0.84 

2.33 

0.73 

1.00 

0.60 

6.1 

7.2 

1 

1.20 

0.84 

1.11 

0.18 

0.24 

0.69 

6.5 

7.7 

2 

1.73 

0.84 

1.66 

0.40 

0.54 

0.72 

6.8 

7.7 

3 

1.62 

0.84 

1.54 

0.34 

0.45 

0.75 

7.6 

7.9 

4 

1.40 

0.84 

1.29 

0.22 

0.29 

0.78 

6.4 

8.0 

5 

1.20 

0.84 

1.14 

0.16 

0.23 

0.82 

— 

8.2 

The  calculated  line  broadening  is  in  excellent  agreement  with  our 
experimental  results.  They  identify  an  elastic  vibrational  dephasing 
component  which  is  relatively  J - independent  and  quite  commensurate  with  our 
estimate.  In  addition  they  find  a  vibrational  state  dependence  to  the 
inelastic  rates  which  we  also  have  identified  in  our  study  of  HD  [P4]  . 
Agreement  at  this  level  of  detail  is  very  encouraging  support  for  the 
empirical  analyses  we  have  applied  to  the  hydrogen  systems,  cf.  Section  III. 
Although  there  have  been  calculations  for  other  H2'.rare-gas  systems,  the 
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current  state  of  our  experimental  characterization  for  these  systems  does  not 
warrant  detailed  comparisons. 

Turning  now  to  self -broadened  H2 ,  we  have  added  significant  new  data  at 
elevated  temperatures.  Unfortunately,  there  are  no  calculations  which  cover 
this  range.  At  room  temperature,  the  ECS  analysis  (discussed  in  Section 
III. B. 2. a)  identifies  an  R-R  dominated,  inelastic  contribution  in  agreement 
with  the  calculation  of  Van  Kranendonk  [R6]  and  a  vibrational  dephasing 
contribution  to  the  linewidth  in  agreement  with  the  calculation  of  Kelley  and 
Bragg  [R17] .  The  model  of  the  latter  study  would  not  predict  a  J  dependence 
for  the  vibrational  dephasing  contribution  to  the  linewidth;  however,  these 
authors  do  predict  an  almost  linear  increase  in  this  component  as  temperature 
is  raised  from  295  to  500  K.  We  note  in  this  regard  that  the  shifting 
coefficient  (which  arises  from  the  same  interaction  as  the  vibrational 
dephasing  portion  of  the  line  broadening)  predicted  by  these  authors  shows  a 
too  rapid  increase  with  T  above  room  temperature,  see  the  discussion  below. 

2)  LINE  SHIFTING 

The  situation  with  respect  to  theoretical  treatments  of  line  shifting  is 
similar  to,  but  even  more  limited  than,  that  described  above.  There  is 
excellent  agreement  between  theory  and  experiment  in  the  case  of  DjiHe,  cf. 
Table  IV. 1.  The  only  other  system  for  which  comparison  can  be  made  is  self- 
broadened  H2 . 

The  work  of  Kelley  and  Bragg  [R17]  attempted  to  account  for  the  overall 
temperature  and  J  dependence  of  the  shifting  coefficients  displayed  in  Fig. 
II.  6.  To  understand  the  context  of  their  analysis,  which  stems  from  the 
earliest  measurements  of  Q-branch  shifting  [R18] ,  we  replot  these  data  versus 
the  fractional  population  of  the  J  state  in  Fig.  IV. 2.  We  have  added  to  this 
plot  the  previously  determined  shifts  at  low  temperatures  [R19] .  The  linear- 
with-population  trend  observed  in  these  data  has  been  attributed  to  resonance 
coupling  [R18].  More  recently  this  notion  has  been  incorporated  into  the 

impact  theory  of  line  broadening  by  Ben  Reuven,  cf.  Ref.  [R20.].  Our  data 

show  this  effect  up  to  1000  K.  The  significant  outliers  (seen  above  room  T) 
to  the  linear-with-population  trends  in  these  data  are  the  Q(0)  line  shift 
values.  What  is  remarkable  in  these  data  is  that  the  coupling  shift 
(determined  from  the  slopes  of  the  straight  lines  in  Fig.  IV. 2)  is 

essentially  temperature  independent  from  80  to  1000  K,  whereas  the  overall 
shift,  as  measured  by  the  zero-population  intercepts  of  the  straight  lines, 
changes  from  strongly  negative  to  strongly  positive  over  this  same 
temperature  range.  Our  current  understanding  of  the  H2:H2  intermolecular 
potential,  in  particular  its  vibrational  state  dependence,  will  not  provide  a 
reasonable  explanation  for  either  of  these  observations.  This  topic  will  be 
discussed  in  some  detail  in  Ref.  [M8].  For  the  present  we  limit  the 

discussion  to  the  observation  that  present  theory  [R17]  differs  significantly 
with  the  overall  shift  coefficient  at  elevated  temperatures,  cf.  Fig.  IV.  3; 
and  thus,  we  have  no  reliable  predictive  basis  for  either  the  J  or  the  T 
dependence  of  the  line  shifting  coefficients. 

3)  SPECTRAL  DISTRIBUTIONS 

Let  us  begin  with  the  observation  that  there  is  no  formal,  theoretical 
derivation  or  justification  for  the  form  of  the  spectral  distribution 
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function  which  we  applied  to  the  pressure  broadened  line  shape  of  the  H2 :Ar 
system.  Our  approach  was  phenomenological  and  derived  from  argument-by¬ 
analogy  from  the  theory  of  Doppler  broadening  and  Dicke  narrowing.  This 
appears  as  an  area  which  would  benefit  from  a  careful  theoretical  analysis. 
It  is  generally  true  that  we  still  lack  a  good  understanding  of  spectral  line 
formation  when  there  are  correlations  among  the  velocity,  speed,  or  phase 
changes  associated  with  collisions.  Fortunately,  the  simple  kinetic  models, 
applied  to  this  point,  do  appear  to  give  reasonable  account  of  the 
observations.  We  lack  a  means  of  predicting  the  "kinetic"  parameters  at 
present,  and  this  is  due  in  part  to  the  lack  of  a  good  formal  theory. 

The  quality  of  our  formal  theory  and  a  demonstration  of  its  utility  for 
first  principle  predictions  are  both  evidenced  by  the  theoretical  work, 
discussed  in-part  above,  on  D2:He.  The  final  subject  of  this  work  [R16]  was 
the  shape  of  the  Q  branch  in  the  density  regime  intermediate  to  the  Doppler 
and  fully  pressure-broadened  limits.  In  this  work,  the  authors  use  recently 
derived  extensions  (which  explicitly  treat  spectral  line  formation)  of  the 
Waldmann- Snider  kinetic  equation.  They  examine  several  approximation 
schemes  for  solution  of  the  generalized  kinetic  equations  and  determine  line 
shape  functions  in  good  agreement  with  our  experimental  results.  As  an 
illustration  of  this  we  show  in  Fig.  IV. 4  a  comparison  between  the  the  soft- 
collision  model  which  we  verified  (see  Section  II. B. 3)  as  appropriate  for  our 
measurements  and  the  numerically  calculated  result  supplied  by  Dr.  Monchick 
[R21].  Again,  we  emphasize  that  the  calculation  is  based  on  an  ab  initio 
potential,  involves  generalized  kinetic  cross-sections  which  appear  in  the 
approximate  solutions  to  the  extensions  of  the  Waldmann- Snider  equation,  and 
employs  S-matrices  calculated  from  close-coupled  quantum  scattering 
calculations.  We  are  encouraged  by  such  calculations  to  suggest  that  a 
generally  useful  basis  for  understanding  and  prediction  spectra  not  only  is 
possible  but  soon  will  be  realized  in  fact. 
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V  THEORETICAL  PREDICTION  OF  LINE  BROADENING  FOR  C0:N2 

The  theoretical/computational  development,  described  in  the  preceeding 
section,  provides  a  powerful  tool  for  the  prediction  of  line  broadening  for 
simple  diatomic  systems.  In  particular,  there  is  a  great  deal  of  interest, 
for  example  on  the  part  of  the  combustion  research  community,  in  the 
prediction  of  Q-branch  spectra  for  the  CO  molecule  in  a  N2  host. 

To  address  this  problem,  the  collision  partner  in  the  previous 
calculation  is  replaced  by  N2  using  the  parameters  contained  in  Table  V.l. 
These  varues  supplement  the  CO  values  given  in  Appendix  A,  Table  2. 


Table  V.l 

N2  and  N2 :C0  Potential  Parameters  for  Line  Broadening 

Calculations® . 


Quantity  Value  Ref. 


fi  (lO'^o  C 

m) 

0.00 

Q  (10'''“  C 

m2) 

-4.67 

[R22] 

(10-22 

J) 

5.03 

[R23] 

(10-“ 

m) 

3.32 

[R23] 

^c-n  <10-^^ 

J) 

5.96 

[R24] 

Oc-n  (10-11 

m) 

3.335 

[R24] 

^o-n  (10-^^ 

J) 

4.72 

[R24] 

‘^c-n  (10-^^ 

m) 

3.235 

[R24] 

In-n  (10-“ 

m) 

1.098 

[R25] 

(10-“ 

J) 

13.0 

[R26] 

<10-“ 

m) 

3.75 

[R26] 

® Intermolecular  potential  parameters  given  here 

are  related  to  the  potential  parameters  of  Refs  [R8-R10]  through 

^nd  jr^^j  +  Ir^^l  = 


On  the  basis  of  our  experience  with  self-broadened  CO  and  N2  [Mil],  we 
estimate  the  inaccuracy  of  the  calculated  line  broadening  to  be  of  the  order 
of  5%.  The  comparison,  displayed  in  Fig.  V.la,  with  the  experimental 
infrared  data  of  Refs.  [R12,  Rll]  shows  an  excellent  level  of  agreement.  The 
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Figure  V.  1 


calculated  values  of  the  Q-branch  line  broadening  as  functions  of  J  and  T 
are  listed  in  Table  V.2. 


Table  V.2 

C0:N2  Q-branch  line  broadening  coefficients  (cm’^/atm)® 


IGoO 


1500  K 


0 

0.07631 

0.04902 

0.03747 

0.02883 

0.02179 

1 

0.07086 

0.04638 

0.03673 

0.02796 

0.02129 

2 

0.0/124 

0.04686 

0.03659 

0.02825 

0.02152 

3 

0.06962 

0.04652 

0.03634 

0.02830 

0.02158 

4 

0.06729 

0.04570 

0.03590 

0.02820 

0.02157 

5 

0.06531 

0.04474 

0.03536 

0.02799 

0.02151 

6 

0.06318 

0.04368 

0.03477 

0.02773 

0.02140 

7 

0.06219 

0.04278 

0.03419 

0.02742 

0.02128 

8 

0.06108 

0.04202 

0.03368 

0.02709 

0.02113 

9 

0.06059 

0.04148 

0.03326 

0.02677 

0.02097 

10 

0.05981 

0.04104 

0.03292 

0.02648 

0.02081 

11 

0.05918 

0.04070 

0.03264 

0.02623 

0.02065 

12 

0  05839 

0.04039 

0.03241 

0.02602 

0.02049 

13 

0.05749 

0.040C'' 

0.03220 

0.02584 

0.02036 

14 

0.05664 

0.03975 

0.03199 

0.02568 

0.02023 

15 

0.05557 

0.03938 

0.03178 

0.02554 

0.02012 

16 

0.05467 

0.03902 

0.03157 

0.02540 

0.02001 

17 

0.05357 

0.03859 

0.03137 

0.02527 

0.01992 

18 

0.05266 

0.03818 

0.03114 

0.02514 

0.01983 

19 

0.05161 

0.03772 

0.03089 

0.02501 

0.01975 

20 

0.05073 

0.03728 

0.03063 

0.02487 

0.01967 

21 

0.04977 

0.03680 

0.03036 

0.02472 

0.01959 

22 

0.04895 

0.03634 

0.03009 

0.02457 

0.01950 

23 

0.04810 

0.03587 

0.02980 

0.02441 

0.01942 

24 

0.04735 

0.03541 

0.02951 

0.02425 

0.01938 

25 

0.04660 

0.03496 

0.02922 

0.02408 

0.01933 

26 

0.04591 

0.03452 

0.02892 

0.02391 

0.01928 

27 

0.04524 

0.03409 

0.02863 

0.02373 

0.01923 

28 

0.04461 

0.03368 

0.02834 

0.02355 

0.01918 

29 

0.04401 

0.03329 

0.02806 

0.02337 

0.01913 

30 

0.04342 

0.03291 

0.02778 

0.02319 

0.01908 

•the  five  significant  figures  represent  only  the  numerical  precision  of  the 
calculation;  the  values  are  considered  accurate  at  the  5%  level 


These  values  are  presented  in  Fig.  V.lb  along  with  a  "MEG"  law  (no 
restriction  on  delta-J)  description  based  on  Eqs .  II.  1  and  II.  2  with  the 
following  parameters: 
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A  =  1.75  (V.l) 
a(To)  =:  0.01204;  (Tq  =  300  K)  (V.2) 
m  =  0.204;  (Tq  =  300  K)  (V.3) 
^(T)  =  (  2  -  0.5131(1000/1)° ■ 3 )  (V.4) 
5(T)  =  {  1  +  0.18(1000/T)°-^M  (V.5) 


The  value  of  A  is  assigned  on  the  basis  of  the  average  of  the  N2  and  CO 
values.  We  use  the  temperature  dependent  forms  for  and  5  in  order  to  give 
a  reasonably  accurate  representation  of  the  calculated  Q-branch  linewidths 
over  the  complete  range  of  J  and  T.  This  calculation  will  be  the  subject  of 
a  forthcoming  publication  [M12]. 
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VI  SENSITIVITY  ANALYSIS 


In  this  section  we  will  summarize  selected  results  from  a  study  [MIA]  of 
the  accuracy  which  can  be  achieved  in  the  determination  of  T  and  P  from  a 
single-shot,  multi-channel  CARS  spectrum.  The  widely  used,  well  studied  N2 
vibrational  (v=0  -»  v=l)  Q-branch  spectrum  is  selected  as  the  basis  for  this 
analysis.  The  approach  is  to  generate  (via  computer  calculations)  a  set  of 
"data"  for  self-broadened  N2  for  values  of  T  and  P  over  the  range  295-1500  K 
and  1-100  atm.  We  select  the  9  combinations  of  295,  1000,  and  1500  K  and  1, 
10,  100  atm  as  our  "known"  set  of  Tg  and  Pg  .  The  spectral  data  are 
determined  from  the  Lest,  currently  available  description  of  the  N2  Q  branch 
[R27,  R28]  ,  i.e.  a  modified  exponential  gap  (MEG)  rate  law  [R29]  for  the 
state-to-state  rates  for  rotationally  inelastic  collisions  combined  with  the 
relaxation  matrix  formalism  for  describing  the  CARS  Q-branch  spectral 
distribution  function.  To  more  closely  model  diagnostic  spectra,  the  CARS 
spectral  distribution  also  is  calculated  for  finite  spectrometer  resolutions 
(0.5  cm' ^  and  2.0  cm' ^  Gaussian  convolutions).  The  three  "data"  sets 
(infinite  resolution,  0.5  cm' ^  and  2.0  cm' ^  convolutions)  are  shown  (as  the 
symbols,  which  may  appear  as  heavy  solid  lines)  in  Figs.  VI. 1-3. 

An  estimate  of  the  uncertainty  limits,  "noise-band",  found  in  the 
measurement  of  a  single-shot,  multi-channel  spectrum  is  derived  from  a 
number  of  literature  sources  [R30-R36] .  Since  a  wide-variety  of  noise 
sources  have  been  considered  in  the  referenced  studies,  we  make  an 
operational  definition  (consistent  with  the  literature)  that  the  frequency 
dependent  noise-band,  AI(u),  is  given  by 

AI(w)  =±7(^2+  (kI(w))2)  (VI. 1) 


with  I(w)  the  CARS  spectral  intensity  in  "counts".  In  this  expression,  we 
ignore  the  shot-noise  (proportional  to  the  yi(a>))  because  it  is  usually  of 
the  order  of  the  noise-proportional-signal  component  only  at  the  l%-of-I„,j; 
level  in  the  spectrum.  We  eliminate  consideration  of  these  low  intensity 
regions  of  the  experimental  spectra  by  setting  an  effective  reading-noise 
limit  o  of 


a 


2%  I„ 


(VI. 2) 


The  noi.se  proportional-to-signal  component  is  set  at  the  x=7.5%  level;  a 
reasonably  typical  value  for  diagnostic  spectroscopic  systems .  This 
representation  of  the  uncertainty  limits  is  appropriate  to  broadband- StoKes 
laser  and  single-mode  Raman-pump  laser  configurations. 

In  order  to  estimate  uncertainty  limits  we  choose  as  our  error 
criterion  the  condition  that  no  point  in  the  "fit"  spectrum  differ  from  the 
"data"  by  an  amount  greater  than  the  error  band  defined  in  Eq.  (VI. 1).  This 
condition  is  illustrated  in  Fig.  VI. A.  We  find  the  uncertainty  limits  by 
searching  P,T- space  around  the  "known"  Pg  and  Tg  to  locate  the  points  at 
which  the  error  condition  is  violated.  At  each  search  point  the  overall 
amplitude  of  the  spectrum  is  adjusted  to  minimize  the  weighted  (Eq.  (VI. 1)) 
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F igure  VI .  1 


mean-square  deviation  of  the  "fit"  relative  to  the  "data".  The  locus  of  such 
points,  the  "error- limit  surface",  for  the  three  data  sets  is  plotted  in 
Figs.  VI. 5-7.  In  these  figures,  the  T,P  limits  are  expressed  as  fractional 
deviations  from  the  reference  Tq,Po. 

Since  this  representation  of  the  error-limit  surface  has  not  previously 
been  presented  in  the  literature,  we  suggest  the  following  interpretation  of 
these  plots.  Consider  an  ensemble  of  single-shot  spectra,  each  with  its  own 
unique  noise  component.  If  we  applied  a  least-mean-squares  fitting  procedure 
to  each  of  these  spectra,  we  would  expect  to  find  a  distribution  of  "best 
fit"  P's  and  T's  in  the  vicinity  of  Pq  and  Tq  .  This  distribution  is  a 
mapping  of  the  spectral  noise  through  the  instrument  response  function  and 
the  CARS  spectral  distribution  function  onto  the  T,P  plane.  Such  a  mapping 
is  shown  schematically  in  Fig.  VI . 8 .  Since  the  most  probable  T  and  P  for  the 
set  of  measurements  should  be  Tg  and  Pg ,  we  would  expect  that  the 
distribution  function  can  be  approximated  by  a  figure  of  revolution  of  a 
statistical  distribution  which  maximizes  at  Tg  and  Pg  .  A  Gaussian 
distribution  is  used  to  generate  Fig.  VI . 8 ,  although  the  actual  distribution 
need  not  be  Gaussian  or  symmetric.  For  the  mapping  displayed  in  Fig.  VI . 8  it 
is  easy  to  show  that  the  most  probable  values  for  constant  P  or  constant  T 
cuts  through  the  surface  are  approximately  as  shown  by  the  solid  and  dash-dot 
lines  on  the  figure.  The  error- limit  surface  then  represents  a  measure  of 
the  limits  of  the  probability  distribution,  e.g.  the  2-sigma  points  are 
suggested  by  the  illustration. 

With  the  interpretation  of  the  proceeding  paragraph  in  mind,  we  have 
performed  least-mean-squares  fits  for  P  at  constant  T  (and  vice  versa)  for 
the  entire,  instrumentally  convolved  "data"  set.  In  each  case  the  amplitude 
is  adjusted  to  minimize  the  weighted  mean-square  deviation  as  discussed 
above.  The  solid  and  dot-dashed  lines  in  Figs.  VI. 6  and  VI . 7  are  the  locus  of 
these  weighted  least-squares  fits  with  P  and  T  held  constant,  respectively. 
We  show  in  Fig.  VI. 9  an  expanded  view  of  the  100  atm,  295  K  panel  from  Fig. 
VI.  6.  It  can  be  seen  from  this  figure  that  the  constant  P  and  constant  T 
curves  approximate  those  we  expected  from  the  statistical  argument  of  the 
proceeding  paragraph.  That  is,  if  we  consider,  as  an  example,  that  subset  of 
fits  enclosed  within  the  limit-surface  in  Fig.  VI . 9  for  which  the  pressure  is 
between  9  and  10%  greater  than  Po=100  atm,  we  will  find  from  this  subset  a 
most  probable  temperature  which  is  2%  lower  than  To=295  K.  If  we  consider 
the  limit-surface  at  1  atm,  1500  K,  shown  expanded  in  Fig.  VI. 10,  we  see  that 
errors  in  P  imply  essentially  no  correlated  error  in  T.  We  also  note  from 
this  last  figure  that  temperature  uncertainties  in  the  range  of  ±5%  would  be 
expected  from  single-shot  CARS  data  on  high  temperature  atmospheric  pressure 
N2  .  This  estimate  is  in  good  agreement  with  the  observations  of  Ref.  [R34] 
on  flames  and  Ref.  [R31]  on  a  high  temperature  furnace. 

We  also  observe,  from  a  comparison  of  Figs.  VI. 5-7,  that  the  effect  of 
the  finite  resolution  of  diagnostic  multi-channel  instruments  is,  by-and- 
large,  simply  to  expand  the  error-limit  surface,  particularly  along  the 
pressure  axis,  but  otherwise  to  preserve  the  basic  shapes  and  orientations  of 
the  limit-surface.  We  note  in  this  regard  that  the  magnitude  of  the  P,T  axis 
in  Fig.  VI . 7  is  approximately  twice  that  in  Figs.  VI . 5  and  VI . 6 . 

An  important  observation  derived  from  the  shapes  of  the  limit-surfaces 
and  locus  of  constant  P  and  constant  T  least-squares  minima  is  that  there  are 
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Figure  VI. 6 
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Figure  VI . 7 
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cases  where  errors  in  P  and  T  are  correlated.  These  cases  are:  1  atm,  295  K; 
10  atm,  1000  and  1500  K;  and  100  atm,  295  K.  It  is  important  to  note  that 
the  exact  form  and  degree  of  correlation  is  a  func  ion  of  the  CARS  spectral 
distributioi.  (i.e.  a  strong  function  of  Pq  and  Tg ) ,  the  spectral  resolution, 
and  the  spectral  distribution  of  the  error  (e.g.  the  weighting  function).  We 
illustrate  sensitivity  to  the  latter  by  the  results  presented  in  Fig.  VI. 11. 
For  the  upper  three  panels  the  weighting  function  includes  a  shot-noise 
component,  y(.5I(to)),  a  noise  floor,  0,3=2  (cts.),  with  (cts.),  and 
the  noise-proportional  to  signal  at  the  7.5%  level  as  previously.  These 
values  are  typical  of  diagnostic  applications  to  laboratory  flames,  cf.  Ref. 
[R34] .  The  lower  three  panels  (in  which  we  indicate  only  the  locations  of 
constant  P  and  constant  T  least  squares  minima)  utilize  a  constant,  unit 
weight,  i.e.  commonly  referred  to  as  "unweighted"  fitting.  Upon  comparing 
this  figure  to  Fig.  VI. 6,  we  see  distinct  changes  in  the  magnitude  and 
orientations  of  the  error-limit  surfaces  and  associated  constant  P  and 
constant  T  least  squares  minima.  In  some  cases  the  degree  of  correlation  in 
the  P  and  T  errors  is  increased,  while  in  others  it  is  decreased. 

We  conclude  this  section  with  a  brief  discussion  of  the  errors  in  T  and 
P  which  arise  if  there  are  systematic  errors  in  our  knowledge  of  the 
parameters  of  the  rate  law  which  is  the  basis  for  describing  the  CARS 
spectrum.  These  (and  other  sources  of  systematic  error)  are  studied  in  Ref. 
[M14] .  We  find,  from  an  analysis  similar  to  the  above,  that  systematic 
errors  in  linewidths  of  the  order  of  5-10%  in  general  lead  to  similar  errors 
in  derived  T  and  P.  We  illustrate  this  type  of  sensitivity  in  Fig.  VI.  12. 
In  this  figure  we  see  that  the  entire  error- limit  surface,  and  the  location 
of  the  most  probable  T  and  P,  are  shifted  of  the  order  of  4-5%  as  a  result  of 
a  (-10%)  variation  in  the  MEG  law  parameter  This  variation  in  the  ^ 
parameter  produces  the  linewidth  variations  shown  in  Fig.  VI. 13;  these 
variations  have  average  magnitudes  of  4.6%  at  295  K  and  3.1%  at  1500  K.  We 
see  from  this  analysis  the  need  for  an  accurate  data  base  to  assure  the 
accuracy  of  spectroscopically  derived  T  and  P. 

The  analysis  above  centered  on  use  of  the  fundamental  (0  -♦  1)  Q  branch 
for  temperature  and  pressure  measurement.  We  believe  that  the  methodology 
developed  in  the  course  of  this  study  provides  a  useful  approach  to 
sensitivity  analyses.  Among  the  many  conclusions  which  we  can  draw  from  this 
this  preliminary  study,  we  consider  the  following  to  be  most  significant: 

-the  nature  of  the  noise  spectrum  and  the  CARS  spectral  distribution 
determines  the  shape  of  the  error  distribution  in  P,T  space  and  thus 
meaningful  answers  to  sensitivity  questions  require  a  specification  of 
both  the  nature  of  the  noise  and  the  spectral  region  from  which  T  and  P 
will  be  determined 

-it  is  very  difficult  to  accurately  determine  pressure  from  the  shape  of 
the  Q  branch 

-accurate  temperature  measurements  for  cases  where  there  is  significant 
uncertainty  in  our  knowledge  of  the  MEG  law  can  not  be  based  solely  on 
the  spectrum  of  the  fundamental  branch. 

We  do  not  regard  any  of  these  observations  3  revolutionary,  rather  we  view 
them  as  more  quantitatively  substantiated  statements  of  our  expectations  as 
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Figure  VI .  12 


Frac.  change  in  P 


to  the  limits  of  ertot.  Further,  we  note  that  the  methodolgy  developed  here 
can  be  used  to  suggest  possible  approaches  to  the  minimization  of  systematic 
errors  in  spectral  measurements  of  T  and  P, 
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VII  CONCLUSIONS 


A  reliable  predictive  basis  has  been  established  for  simple  molecular 
systems  such  as  self -broadened  CO  and  N2 .  The  very  detailed  studies  made  on 
the  CO  system,  in  particular,  support  the  use  of  a  relaxation  matrix,  based 
on  rates  for  state-to-state  rotational  energy  transfer,  for  describing  the  0- 
branch  spectrum.  These  studies  also  illustrate  the  limits  of  our  ability  to 
specify  the  relaxation  matrix.  We  observe  that  as  we  attempt  to  account  more 
accurately  for  the  elements  of  the  matrix,  we  encounter  a  more  complex  and 
(currently)  less  predictable  J  and  T  dependence  in  these  elements. 
Fortunately,  the  MEG  law,  in  spite  of  inaccuracies  with  respect  to  a  number 
of  significant  details,  does  give  a  simple  and  remarkably  useful 
representation  of  the  collapse  of  the  Q  branch  as  a  function  of  P  and  T. 

Ultimately,  the  validity  of  the  predictive  model  must  be  judged  in  the 
context  of  our  application,  i.e.  considering  the  diagnostic  approach,  the 
accuracy  required,  ; nd  the  ranges  of  extrapolation  in  pressure,  temperature, 
and  species  concentration  over  which  we  will  try  to  apply  the  predictive 
basis.  It  is  the  fact  that  we  almost  always  must  extrapolate  our  predictive 
equations  into  untested  ranges  which  requires  that  we  know  the  line-formation 
mechanisms  and  their  basic  dependence  on  the  "system"  variables. 

It  is  clear  from  our  studies  that  the  reliability  of  tne  predictive 
basis  would  be  greatly  enhanced  by  first  principle  calculations  of  the  off- 
diagonal  elements  of  the  relaxation  matrix.  Since  these  are  not  strictly 
state-to-state  rates  this  will  require  that  line -broadening- type 
calculations  be  extended  to  reveal  the  "distribution  of  relaxation"  as  a 
function  of  the  intermolecular  in. tractions  and  the  magnitudes  of  angular 
momentum  and  energy  transfer. 

The  need  for  a  first  principle  understanding  increases  significantly 
when  we  include  the  composition  variable  in  our  systems.  This  situation 
arises  Vecause  it  is  very  difficult  to  empirically  study  the  spectrum  of  the 
diagnostically  interesting,  minority,  molecular  species  in  the  presence  of 
the  foreign  gas  host.  The  difficulty  involves  not  only  making  accurate 
spectroscopic  measurements  on  a  minority  species,  but  also  problems 
associated  with  producing  stable,  accurately  known  reference  samples  having 
the  required  ranges  of  concentration,  pressure,  and  temperature.  In 
addition,  no  matter  what  limitations  are  set  by  the  experimental  difficulty 
of  an  empirical  validation-study,  diagnostic  applications  will  (by 
definition)  almost  always  involve  significant  extrapolations  of  some  or  all 
of  the  system  variables.  Improved  first  principle  understanding  would  be 
extremely  valuable. 

For  molecules  such  as  N2  and  CO,  the  development  of  a  predictive  basis 
for  the  mixed-gas  systems  also  will  require  significant  improvements  in  the 
sensitivity  and  accuracy  of  the  spectroscopic  measurement  techniques  applied 
in  model-validation  studies.  High  resolution  CARS  spectroscopy  appears  to  be 
a  reasonable  approach  to  provide  some  of  this  capability.  Additionally,  the 
fact  that  we  now  have  accurately  known  susceptibilities  for  certain  systems 
can  provide  a  means  of  discerning  systematic  errors  in  multi-channel 
(diagnostic)  instruments.  The  last-mentioned  instruments,  so  calibrated, 
then  can  be  used  on  reference  samples  to  provide  an  important  level  of 
validation . 
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The  spectrum  of  the  hydrogens  is  very  appealing  from  the  diagnostic 
point  of  view,  in  that  the  resolved,  intense  Q  branch  affords  excellent 
measures  of  temperature  and  pressure.  We  note  also  that  the  large  pressure- 
shift  to  pressure -width  ratio  characteristic  of  the  hydrogens,  especially  in 
foreign  gas  hosts,  suggests  a  unique  approach  to  pressure  measurements  which 
overcomes  the  large,  correlated  ambiguities  in  pressure  and  temperature 
measurements  which  can  be  encountered  when  using  the  Q  branch  of  molecules 
such  as  CO  and  N2 . 

On  the  basis  of  our  studies  of  the  hydrogens,  we  see  that  a  reasonable 
conceptual  basis  for  description  of  these  systems  is  emerging  and  that 
empirical,  predictive  formulae  can  be  used  for  some  important  parameters. 
However,  there  remain  fundamental  questions  with  regard  to  the  physics  of 
line  formation  as  well  as  problems  with  respect  to  reliable  prediction  of 
spectroscopic  parameters.  The  hydrogen  diatomics,  although  quite  tractable 
in  terms  of  the  spectroscopy,  require  significant  levels  of  theoretical 
development,  e.g.  new  kinetic  equations,  improved  intermolecular  potentials 
at  short  range  including  vibrational  state  dependence,  and  computationally 
manageable  quantal  scattering  methods.  Producing  requisite  reference  samples 
also  can  pose  significant  problems.  The  study  of  the  hydrogens  affords  many 
opportunities  for  fundamental  advances  in  science  and  has  significant  promise 
for  utility  in  diagnostic  application. 

We  note  that  an  element  common  to  all  fundamental  and  model-validation 
studies  is  the  necessity  of  producing  the  required  "reference  samples". 
These  must  have  accurately  known  conditions  of  temperature,  pressure,  and 
species  concentration,  selected  either  to  be  reasonably  representative  of  the 
diagnostic  environments  of  interest  or  to  critically  test  an  element  of  the 
predictive  theory.  In  this  regard,  we  suggest  that  shock  tubes  may  serve  as 
sources  of  these  reference  samples.  Use  of  these  sources  will  demand 
accurately  calibrated  multi-channel  spectroscopy;  which,  as  we  indicated 
above,  now  appears  within  reach  for  both  low  and  high  resolution 
applications.  The  combination  of  a  shock  tube  source  with  a  "reference 
quality"  spectroscopic  instrument  could  provide  important  checks  on  possible 
systematic  errors  in  pressure,  temperature,  and  species  concentration  for  a 
wide-range  of  diagnostic  problems. 

Finally,  we  observe  that  the  accuracy  required  in  the  data  base  in  order 
to  avoid  significant  systematic  errors  should  be  assessed  in  terms  of  some 
form  sensitivity  analysis.  Such  analyses  will  guide  the  selection  of  the 
range  of  conditions  under  which  validation  experiments  are  required. 
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Abstract 


We  test  the  ability  of  the  Energy  Corrected  Sudden  (ECS) ,  Modified 
Exponential  Gap  (MEG)  and  the  Statistical  Power-Exponential  Gap  (SPEC)  rate 
laws  to  characterize  line  broadening  and  line  interference  in  the  CO  Q-branch 
at  295  K.  All  three  rate  laws  fit  the  experimental  linewidth  data.  The  ECS 
law  is  found  to  predict  too  much  spectral  collapse.  The  MEG  and  SPEC  laws 
both  adequately  model  spectral  collapse,  but  with  different  implications 
about  the  role  of  dipolar  and  quadrupolar  symmetry  forces  in  CO: CO  line 
broadening.  From  semiclassical  calculations  of  CO  linewidths,  we  conclude 
that  the  SPEC  law  with  a  restriction  to  even  AJ  changes  is  the  more 
physically  correct  model. 


I .  Introduction 


Numerous  studies  of  the  pressure  and  temperature  dependence  of  Q- 
branch  spectra  of  simple  molecular  systems  [1]  show  that  collisional 
broadening  and  spectral  line  interference  in  the  Q-branch  can  be  described 
using  the  relaxation  matrix  formalism.  [2-5],  When  vibrational  dephasing 
(for  IR  and  Raman  studies)  and  elastic  reorientation  (for  IR  studies)  can  be 
neglected,  the  elements  of  the  relaxation  matrix  are  expressible  in  terms  of 
state-to-state  rates  for  rotational  energy  transfer.  This  has  focused  the 
attention  of  many  studies  on  tests  of  the  ability  of  various  rotationaliy 
inelastic  rate  law  models  [6],  which  provide  a  simple  representation  of  the 
relaxation  matrix,  to  yield  an  accurate  characterization  of  the  observed 
spectra. 

The  present  study  is  motivated  by  the  desire  to  test  the  ability  of 
three,  widely  applied,  rate  law  models  to  predict  the  self -broadened  CO  Q- 
branch  spectrum.  Aside  from  the  potential  implications  for  diagnostics 
applications,  line  interference  in  the  CO  Q  branch  is  of  particular  interest 
for  a  number  of  reasons.  First,  spectroscopic  evidence  [7-12]  gives  a  clear 
indication  that  CO  line  broadening  shows  little  vibrational  state 
dependence,  thus  supporting  the  use  of  a  rate  law  model  to  parameterize  the 
relaxation  matrix.  Second,  neighboring  transitions  in  the  CO  Q  branch  when 
overlapped  may  interfere,  thereby  causing  the  spectrum  to  collapse.  This 
degree  of  line  interference  should  provide  a  good  test  of  rate  law  models. 
The  situation  for  CO  contrasts  with  N2  for  which  the  line  interference  is 
obscured  as  the  observed  spectrum  is  the  superposition  of  two  independent  Q 
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branches  from  the  ortho  and  para  molecular  species.  Finally,  there  have  been 
inelastic  rate  law  ana}  ,'ses  of  CO  linewidth  data  [11-13];  however,  no  rate 


law  model  has  been  adequately  tested  against  the  observed  spectral  collapse. 

In  this  paper,  we  present  experimental  results  from  a  series  of 
studies  on  the  pressure  dependence  of  the  CO  Raman  Q-branch  from  0.5  to  6.0 
atm  at  room  temperature.  From  these  studies  we  report  a  J- independent  shift 
coefficient,  improved  line  broadening  coefficients  for  the  lower  J 
transitions  (J=0-18) ,  and  accurate  line  broadening  coefficients  for  the 
higher  J  transitions,  up  to  J=24. 

We  found  from  our  analyses  that  linewldths  and  collapsed  spectra  were 
both  adequately  modelled  by  either  the  Modj.^.  i  Exponential  Gap  (MEG)  [14]  or 
Statistical  Polynomial -Exponential  Gap  (SPEC)  [6,15]  rate  laws,  but  with 
different  implications  as  to  the  relative  propensity  for  dipolar 
(unrestricted  AJ)  and  quadrupolar  (even  AJ  only)  transitions.  The  Energy 
Corrected  Sudden  (ECS)  [6,16]  scaling  law,  previously  reported,  can  be 
discounted  as  it  predicts  too  much  spectral  line  interference.  From  the 
results  of  a  semiclassical  calculation  of  CO  line  broadening,  we  conclude 
that  the  SPEC  rate  law  model,  with  a  restriction  to  quadrupolar  selection 
rules  governing  rotational  state  changes,  provides  the  more  physical 
description  of  the  CO  Q-branch  spectrum  at  295  K. 

II .  Experimental  Results 

The  data  reported  here  are  derived  from  three  experimental  studies  of 
the  self-broadened  Q  branch  of  CO  at  room  temperature.  The  earliest  of  these 
[10]  employed  inverse  Raman  spectroscopy  with  cw-lasers  and  rletermined  line 
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broadening  out  to  J=18  and  a  line  shifting  coefficient  for  Q(14).  A  second 
series  of  experiments  [11],  conducted  with  a  pulsed  Raman  pump  laser  in  a 
quasi -cw  inverse  Raman  spectroscopic  system,  extended  line  broadening  out  to 
J=24  and  measured  shift  coefficients  for  Q(J)=18-22.  A  final  set  of 
measurements,  using  the  latter  experimental  system,  were  made  to  accurately 
determine  the  spectrum  of  the  partially  collapsed  Q  branch  over  a  large 
frequency  range  (above  and  below  the  bandhead)  at  2.9  and  6.0  atm.  These 
experimental  systems,  the  associated  spectral  resolution  and  frequency 
accuracy,  and  the  sample  cells  and  optical  configurations  have  been 
described  in  previous  publications  [10-12]. 

A.  Line  broadening  coefficients 

Using  the  fitting  procedure  discussed  in  Refs.  [11]  and  [12],  we 
determined  the  line  broadening  coefficients  given  in  Table  1  and  displayed 
in  Fig.  1.  Error  estimates,  indicated  in  parentheses,  are  typically  on  the 
order  of  2-5%.  The  low  J  transitions  (0-3)  are  strongly  blended  leading  to  a 
larger  uncertainty.  The  high  J  transitions,  due  to  their  relative  weakness, 
also  have  larger  uncertainties.  It  should  be  emphasized  that  the  values 
extracted  through  our  fitting  procedure  are  independent  of  any  rate  law 
model . 

B.  Line  shifts 

To  accurately  predict  the  collapse  of  a  Q  branch  it  is  important  to 
have  a  reliable  value  of  the  pressure  shift  coefficient  [11,17].  Since 
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collisional  collapse  has  che  effect  of  "pulling”  Q-branch  lines  toward  the 
most  probable  J,  it  is  important  that  line  mixing  be  small  if  one  is 
interested  in  accurate  frequency  shift  determination  [11].  The  measured 
transitions  J=18-22  display  a  sufficiently  high  signal-to-noise  ratio,  and 
have  small  enough  line  mixing  for  the  pressure  range  used  (0.5  -  2.9  atm), 
that  accurate  values  of  the  line  shifting  coefficient  can  be  determined. 
For  the  range  J=18-22,  there  is  no  discernible  J  dependence  of  the  pressure 
shift.  We  note,  that  this  observation  is  consistent  with  measurements  [17] 
and  calculations  [18]  for  self  broadened  N2  ,  and  with  the  results  of  our 
calculation,  reported  below.  A  weighted  linear  regression  of  the  shifts 
versus  pressure  yields  a  value  of  the  CO  Q-branch  shift  coefficient,  -0.0032 
±  0.0003  cm" ^ /atm  at  298  K. 

III.  Rate  Law  Analysis  of  CO  0- Branch  Data 

The  rate  law  models  we  investigated  are  the  Energy  Corrected  Sudden 
(ECS)  [16],  the  Modified  Exponential  Gap  (MEG)  [14],  and  the  Statistical 
Power-Exponential  Gap  (SPEG)  laws  [15],  All  three  of  these  rate  law  models 
have  been  applied  in  the  analysis  of  partially  or  fully  collapsed  Q-branch 
spectra  with  varying  degrees  of  success. 

The  ECS  scaling  law  was  used  by  Bel  Bruno  et  al.  [13]  to  derive 
inelastic  state-to-state  rates  from  CO  infrared  linewidths.  The  ECS  scaling 
law  also  has  been  applied  to  line  broadening  and  collapse  of  the  isotropic 
Raman  Q-branch  spectrum  of  N2  with  moderate  success  [18]. 

The  MEG  law,  which  has  been  the  most  successful  in  characterizing  the 
isotropic  Raman  Q-branch  spectrum  of  N2  [19],  has  been  applied  to  the 
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parameterization  of  the  J  and  temperature  dependence  of  the  CO  Q-branch 
linewidths  [12].  Preliminary  analyses  of  CO  Q-branch  collapse  using  this 
law  have  been  reported  [11,12]. 

The  SPEC  law  has  been  employed  successfully  in  the  description  of  the 
infrared  3  ^-2  band  of  pure  N2O  [20]  and  to  the  2/2  band  of  CO2  ,  both 
self-  [21]  and  N2 -broadened  [22],  It  has  also  been  applied  in  the 
description  of  the  NO  isotropic  Raman  Q-branch  spectrum  [23] .  No  analyses  of 
CO  linewidths  or  spectra  using  the  SPEC  law  have  been  reported. 

A.  The  ECS  Scaling  Law  Analysis 

We  turn  now  to  predictions  of  the  Raman  isotropic  Q  branch  based  upon 
an  ECS  scaling  law  analysis  of  the  CO  infrared  linewidths  [13].  According 
to  this  analysis  two  raechanisms  termed  R-T  R-R  processes,  contribute  to 
inelastic  rates  in  CO.  The  former  is  essentially  a  rotation  to  translation 
energy  exchange  process  for  molecule-"effective”  atom  interactions  and  the 
latter  a  rotational-quantum  number  perserving  exchange  process.  The  R-T 
transfer  rates  are  determined  by  a  set  of  fundamental  base  rates  for  J-+0 
transitions  and  an  associated  scaling  length.  The  fundamental  base  rates  for 
the  R-T  transfer  process  are  expressed  in  terms  of  an  exponential  gap  law 
with  two  parameters.  The  R-R  transfer  rates  are  determined  by  a  single  AJ=1 
exchange  rate  and  an  associated  scaling  length.  We  use  the  parameters  given 
in  Ref,  [13]  along  with  the  definitions  and  equations  given  in  Refs.  [16,24- 
26]  to  calculate  the  ECS  rates. 

We  plot  in  Fig.  1,  along  with  the  experimental  broadening 
coefficients,  the  total  inelastic  rates  (from  both  R-T  and  R-R  processes)  as 
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given  by  the  ECS  analysis.  It  can  be  seen  that  the  ECS  law  gives,  overall, 
reasonably  good  agreement  with  measured  linewidths.  The  ECS  law,  when  used 
to  predict  the  collapsed  Q-branch  spectra  (a  comparison  with  experimental 
results  at  3.0  atm  is  shown  in  Fig.  2),  is  found  to  overestimate  the  degree 
of  collapse.  This  represents  an  example  of  the  fact  that  the  ability  of  a 
rate  law  to  model  linewidths  does  not  necessarily  indicate  its  ability  to 
predict  spectral  collapse. 

Although  the  ECS  analysis  of  Bel  Bruno  et  al.  did  not  discriminate 
between  dipolar  and  quadrupolar  selection  rules,  we  have  subsequently 
performed  an  ECS  analysis  of  CO  linewidths  based  on  quadrupole  selection 
rules,  using  either  an  exponential  gap  or  a  polynomial  gap  law  for  the  base 
R-T  rates  (no  R-R  rates)  .  We  again  find  that  the  ECS  law  significantly 
overestimates  the  degree  of  collapse  and  thus  it  appears  inadequate  for 
describing  the  CO  relaxation  matrix.  The  ECS  ].aw  predicts  inelastic  rates 
near  the  diagonal  which  are  too  large  and  too  little  a  propensity  for  multi¬ 
quanta  collision  processes. 


B.  The  MEG  Law 


The  MEG  law  used  in  Refs.  [11]  and  [12]  is  an  exponential  gap  law 
with  an  energy  correction  term  applied  to  the  initial  level.  We  have 
generalized  the  MEG  law  by  applying  a  similar  energy  correction  to  the  final 
state  following  the  prescription  of  Koszykowski  et  al.  [14].  This  MEG  law 
expresses  the  rate  for  upward  rotational  transitions  from  state  J  to  state 
J '  as 


Wj  .  j  =  P  a 


fl  +  A  E(J')/kT5'  ]2 

1  +  A  E(J)/kT5 

[l  +  A  E(J')/kT  J 

1  +  A  E(J)/KT 
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exp(-  (.) 

with  a,  fi,  S',  and  6  being  adjustable  parameters,  A  a  constant  (a  measure  of 
the  duration  of  a  collision),  E(J)  tb-’  energy  of  the  rotational  level  J,  P 
pressure,  T  temperature,  and  k  Boltzmann's  constant.  We  note  that  in  the 
limit  that  6'=  1,  we  recover  the  form  of  the  MEG  law  given  in  [12,14], 
Previously  the  value  A  =  2  was  assigned  based  upon  estimates  of  effective 
interaction  length  and  distance  of  closest  approach  for  CO:CO  collisions  at 
295  K  [12],  We  note,  however,  that  the  value  of  A  is  not  critical  in  the 
ability  of  the  MEG  law  to  fit  the  data. 

A  best  fit  of  the  MEG  law  to  all  the  experimental  data  (both 
linewidths  and  partially  collapsed  spectra)  gives  a  =  0.01395  (em'^/atm),  = 
1.351,  S'  =  6  =  1.17,  without  any  distinction  between  even  and  odd  AJ 
transitions.  In  Fig.  1  the  MEG  law  is  seen  to  reproduce  the  linewidths  at  a 
level  that  it  consistent  with  experimental  error  at  all  J .  The  results 
obtained  with  this  form  of  the  MEG  law  and  6  =  S'  yield  only  slightly  better 
fits  to  linewidths  and  predictions  of  collapse  than  the  MEG  law  discussed 
previous!’'  in  Refs.  [11,12]  with  S'  =  1. 

In  Fig.  3  a  and  b  we  compare  MEG  law  predictions  with  experimentally 
measured  spectra  at  3.0  and  6.0  atm.  The  residuals,  which  we  assess  to  be 
within  experimental  reproducibility,  demonstrate  that  the  MEG  law  yields  a.n 
excellent  characterization  of  the  Q-branch  collapse.  It  has  been  found  that 
the  degree  of  collapse  observed  for  the  CO  spectrum  is  poorly  predicted  by  a 
MEG  law  with  the  restriction  to  even  rotational  quantum  number  changes, 
which  might  be  expected  for  an  "almost  homonuclear"  diatomic. 
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From  the  MEG 


law  analysis,  we  would  be  led  to  the  conclusion  that  there  is  roughly  equal 
propensity  for  dipolar  and  quadrupolar  transitions  in  CO: CO  collisions. 

C,  The  SPEC  Law 


The  SPEC  law  is  a  statistical  rate  law  model  which  has  its  origins  in 
surprisai  analysis  [6,15] .  This  law  expresses  the  rate  for  an  upward 
rotational  transition  from  state  J  to  state  J'  as 


'j '  J 


=  P  a 


E(J')  -  E(J) 


exp 


-  E(J))^ 
kT 


(2) 


with  a,  /9,and  6  adjustable  parameters  and  By  is  the  ground  state  rotational 
constant,  the  remaining  symbols  are  as  defined  in  section  IIB  above. 

It  was  found  that  the  overall  best  fit  to  the  experimental  data,  both 
linewidths  and  collapsed  spectra,  is  given  by  a  =  0.06428  (cm'^/atm),  S  = 
0.2646  and  p  =  1.295  with  the  collisional  selection  rule  AJ  =  =  even. 
The  SPEC  law  is  also  found  to  fit  the  experimental  linewidths  totally 
consistent  with  experimental  error,  see  Fig.  1. 

In  Figs.  4  a  and  b  we  shew  the  excellent  agreement  of  the  SPEC  law 
predictions  with  our  experimental  spectra  at  3.0  and  6.0  atm.  A  SPEC  law 
without  the  restriction  to  even  AJ  is  found  to  significantly  overestimate  the 
degree  of  collapse  of  the  CO  Q  branch.  From  this  analysis,  we  would  conclude 
that  rotational  state  changes  in  C0:C0  collisions  are  uorainated  by 
quadrupolar  selection  rules.  This  is  in  contrast  to  conclusions  based  on  the 
MEG  law  analysis.  Without  further  information  it  was  not  possible  to  choose 
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between  the  MEG  and  SPEC  laws  and  their  different  implications  about 
rotational  state  changes  in  CO: CO  collisions. 

IV .  Semiclassical  Linebroadening  Calculations 

To  discern,  on  a  physical  basis,  between  the  SPEC  and  MEG  rate  law 
models,  we  have  turned  to  theory  with  the  intention  to  determine  the 
relative  importance  of  dipolar  and  quadrupolar  forces  in  CO; CO  collisions. 
We  use  the  Robert  and  Bonamy  (27-33]  formulation  of  molecular  line 
broadening  to  calculate  CO  self-broadening  coefficients.  In  our 
calculations,  all  terms  from  a  fourth-order  expansion  of  the  atom-atom 
potential  are  explicitly  included  [27-29,37],  The  CO  intermolecular 
interaction  parameters,  taken  from  the  literature,  are  given  in  Table  2. 
Full  details  of  these  calculations  will  be  presented  elsewhere  [43] . 

Our  calculations  show  excellent  agreement  when  compared  with 
experimental  infrared  R-branch  [7,44]  and  Raman  Q-branch  [12]  data  over  a 
large  range  of  J  (0-30)  and  T  (100  -  1500  K)  for  both  infrared  and  Raman  Q- 
branch  data.  Agreement  over  this  range  of  J  and  T  indicates  that  this  CO: CO 
potential  is  very  reliable  for  line  broadening  calculations. 

To  determine  the  role  of  the  quadrupolar  and  dip-^lar  symmetry  forces 
in  CO; CO  collisions,  we  repeated  »-he  calculations  omitting  terms  arising  from 
dipolar  symmetry  forces.  We  find  that  the  dipolar  forces  make  a  small 
contribution  (4-8%)  to  Q-branch  linewidths  for  the  low  and  intermediate  J 
(J=0-14)  transitions  and  have  a  larger  contribution  for  high  J  transitions 
(J>13)  ,  with  =  18%  of  the  Q(24)  transition  lir.  dth  being  attributable  to 
these  forces.  Because  of  the  large  interline  spacing,  the  latter  region  of 
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the  Q  branch  exhibits  relatively  little  line  interference;  in  contrast,  at 
low  J  where  collapse  is  most  important,  we  infer  from  these  calculations  that 


the  CO  state-to-state  transitions  are  dominated  by  quadrupolar  selection 
rules.  Finally,  we  note  that  this  calculation  provides  only  total  inelastic 
collision  cross  rates,  not  individual  state-to-state  rates;  so  it  cannot  be 
used  to  directly  determine  the  off-diagonal  elements  of  the  relaxation 
matrix . 

V .  Conclusion 


The  comparision  of  the  MEG  and  SPEC  law  for  modeling  of  the  CO  Raman 
Q-branch  represents  an  important  example  of  the  fact  that  it  is  difficult  to 
draw  reliable  conclusions  about  collision  dynamics  from  rate  law  fits  to 
linewidths  and  spectra  alone.  To  illustrate  the  difference  between  these  two 
rate  law  models,  we  plot  in  Fig.  5a  and  5b,  respectively,  rotationally 
inelastic  state-to-state  rates  based  upon  the  best  MEG  and  the  SPEG  laws.  As 
can  be  seen  from  this  figure,  the  distribution  of  state-to-state  rates 
given  by  the  two  rate  laws  is  very  different.  Only  the  SPEG  law  analysis  is 
consistent  with  the  dominant  role  of  quadrupolar  symmetry  forces,  as  found 
from  our  calculations.  This  observation  is  also  in  accord  with  experimental 
[45]  and  theoretical  [46,47]  studies  of  CO;He  [47]  and  other  "almost 
homonuclear"  diatomics  in  collision  with  rare  ga'-'-  45].  Our  conclusions 
could  be  tested  by  measuring  the  relative  nj=^  and  AJ=2  state-to-state 
rates.  Finally,  we  note  that  this  study  point<^  '^''t  the  need  for  detailed 
information  about  state-to-state  rates  for  physically  reliable  modeling  of 
high  pressure  Q-branch  spectra. 
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Figures 

1.  Isotropic  Raman  Q-branch  pressure  broadening  coefficients  as  a 
function  of  rotational  quantum  number  J.  The  filled  circles, 
with  error  bars,  are  the  experimental  values  from  Table  1, 

The  dotted  line  is  the  ECS  prediction.  The  dashed  line  is  the 
SPEC  law  fit  given  in  section  III-C,  The  solid  line  is  the  MEG 
law  fit  given  in  section  III-B.  (see  text) 

2.  The  ECS  scaling  law  prediction  of  the  CO  Q-branch  spectrum  at  3.0  atm  and 
295  K  (solid  line)  d  to  experimental  spectra  (dots) .  The 

residuals  (observed  minus  calculated)  are  offset  for  visibility. 

3.  The  MEG  law  prediciton  of  the  CO  Q-branch  spectrum  (solid  line)  compared 
with  experimental  spectra  (dots)  at  (a)  3.0  and  (b)  6.0  atm.  The 
residuals  (observed  minus  calculated)  are  offset  for  visibility. 

4.  The  SPEC  law  prediciton  of  the  CO  Q-branch  spectrum  (sc  id  line)  compared 
with  experimental  spectra  (dots)  at  (a)  3.0  and  (b)  6.0  atm.  The 
residuals  (observed  minus  calculated)  are  offset  for  visibility. 

5.  Rotationally  inelastic  state-to-state  i.ates  as  a  function  of  intial  and 
final  quantum  number  as  given  by  (a)  the  best  MEG  law  model  (5 '=6, 
unrestricted  AJ) ,  and  (b)  the  best  SPEG  law  model  (AJ  =  even  only) . 
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Table  1 


Experimenatal  values  of  CO  self -broadening  coefficients  at  295  K 
(HWHM  in  cm" ^ /atm,  (uncertainty  in  last  digit)) 


J 

0 

0.086(10) 

1 

0.0789(50) 

2 

0.0744(30) 

3 

0.0712(30) 

4 

0.0690(20) 

5 

0.0672(10) 

6 

0.0659(10) 

7 

0.0650(10) 

8 

0.0640(10) 

9 

0.0634(10) 

10 

0.0620(10) 

11 

0.0615(10) 

12 

0.0599(10) 

13 

0.0593(10) 

14 

0.0589(10) 

15 

0.0574(10) 

16 

0.0561(10) 

17 

0.0560(10) 

18 

0.0549(10) 

19 

0.0543(10) 

20 

0.0538(10) 

21 

0.0526(10) 

22 

0.0511(10) 

23 

0.0510(40) 

24 

0.0472(80) 
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Table  2 


Intermolecular  Parameters  for  Line  Broadening 
Calculations® . 


Quantity 

Value 

Ref . 

H  (lO'^o  C 

m) 

.367 

[38] 

Q  Cm) 

-6.00 

[38] 

^c-o  (10-"=^ 

J) 

8.50 

[39] 

°o-c  (10'“ 

m) 

2.95 

[39] 

e,.c  (10-22 

J) 

7.07 

[40] 

^c-c  (10-^' 

m) 

3.35 

[40] 

(10-^^ 

J) 

1  .lb 

[41] 

(10-'* 

n) 

3.15 

[41] 

ic-o  (10-“ 

m) 

1.13 

[36] 

(10-^^ 

J) 

13.8 

[42] 

^iso  (10-“ 

m) 

3.7b 

[42] 

“Intermolecular  potential  parameters  given  here  are  related  t' 
the  potential  parameters  of  Refs  [27-30]  through 

/3 

and  |ri, |  +  [r^^ |  = 
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Abstract 

We  report  observations  of  inhomogeneous  broadening  in  the  vibrational  line 
profiles  of  a  gas  in  the  "impact"  density  regime.  In  measurements  up  to  27  amagat 
(where  the  spectra  are  dominated  by  collision  broadening),  non-Lorentzian,  asym¬ 
metric  features  are  observed  in  Raman  {2-branch  transitions  of  H2  dilute  in  a  heavy 
pertiu’ber  gas.  We  compare  these  measurements  with  an  inhomogeneous  line-profile 
model  based  on  collisional  shifts  that  strongly  depend  on  the  (thermal)  speed  of  the 
H2  radiator.  Quantitative  agreement  is  obtained  only  when  speed-changing  collisions 
are  considered,  which  result  in  spectral  line  narrowing. 
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The  collisional  broadening  of  an  isolated  rotation-vibration  transition  of  a  gas 
is  considered  to  be  in  the  impact  regime^  when  (a)  the  collision  duration  is  small 
compared  to  the  inverse  of  the  collisional  width,  (b)  the  latter  is  large  compared  to  the 
Doppler  width,  and  (c)  the  collisional  width  and  shift  have  no  significant  dependence 
on  relative  molecular  speeds.  Such  conditions  have  been  shown  theoretically^  to 
produce  Lorentz  profiles.  Experimental  studies^  of  collisional  broadening  in  pure  H2 
have  successfully  employed  these  assumptions  for  densities  well  above  the  Dicke 
minimum'*  for  the  Doppler  effect  (-2.5  amagat  for  pure  H2)  and  below  densities  where 
the  finite  duration  of  collisions  becomes  important  (-100  amagat).  However, 
measurements  of  the  temperature  dependence  of  H2  collisional  shift  coefficients^-^ 
indicate  that  such  shifts  depend  significantly  on  relative  speed,  in  apparent  violation  of 
(c).  In  considering  spectra  at  densities  in  the  Doppler  regime,  Berman, ^  Ward  et  al.P 
and  Nienhuis^  demonstrated  that  such  speed-dependent  shifts  (and/or  speed- 
dependent  collision  widths)  can  lead  to  asymmetries  and/or  additional  inhomo¬ 
geneities  in  the  line  profile.  Berman,  Pickett^  and  Luijendijk^O  showed  theoretically 
that  these  effects  generally  persist  to  higher  densities,  i.e.,  to  the  collision-broadened 
regime.  While  examples  of  speed-dependent  line  profiles  in  the  Doppler  regime  have 
been  reported  in  atomical  Ca  and  rovibrational  HF  and  HCl  spectra, to  our 
knowledge  none  have  been  observed  in  a  purely  collision-broadened  line. 

In  this  Letter,  we  report  observations  of  line  profiles  in  Raman  2-branch 
transitions  of  H2  perturbed  by  Ar  that  exhibit  asymmetries  and  Gaussian-like  cores  at 
densities  normally  considered  to  be  in  the  impact  regime.  For  dens'  f ■om  3.6  to  14 

amagat  at  room  temperature  the  non-Loren tzian  features  do  not  vary,  t  iiough  the 
width  of  the  profile  increases  linearly.  Spectra  measured  at  constant  density  (-14 
amagat)  show  a  marked  increase  in  the  non-Lorentzian  features  with  increasing 
temperature.  The  features  occur  at  densities  well  above  the  Dicke  minimum 
(-1  amagat  for  H2-Ar)  and  yet  below  those  for  which  finite  collision-duration 
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effects *3  or  collision-induced  scattering*^  contribute  to  the  Raman  spectrum.  We 
show  that  the  anomalies  can  be  explained  by  inhomogeneous  broadening  that  arises 
from  strong  variations  in  the  collisional  shift  with  radiator  speed.  This  interpretation 
is  supported  by  a  line-profile  model  based  on  speed-dependent  shift  cross-sections 
derived  from  measurements  of  the  Q-branch  line  shifts  of  the  H2-H2  and  H2-Ar 
systems  at  different  temperatures.  The  model  also  incorporates  the  effects  of  speed¬ 
changing  collisions,  which  coherendy  narrow  the  (speed-dependent)  inhomogeneous 
part  of  the  profiles.  Such  collisions  are  most  likely  when  the  masses  of  radiator  and 
perturber  molecules  are  similar,  explaining  why  these  inhomogeneities  are  less 
important  in  self-broadened  gases. 

T^e  high -resolution  inverse  Raman  (IRS)  experiment* ^  employs  a  pulse- 
amplified,  tunable  pump  laser  (~22  ns  pulse  duration)  and  a  cw  probe  laser,  providing 
a  spectral  resolution  of  45  MHz  (FWHM).  Q-branch  transition  frequencies  were 
obtained  using  a  travelling  comer-cube  wavemeter.  The  highest  peak  intensity 
produced  at  the  focus  of  the  beams  was  kept  below  ~2  GW/cm^  to  eliminate  Stark 
effects,  which  have  been  known  to  produce  asymmetries.  *6  The  sample  gases  were 
contained  in  high-pressure,  high-temperature  furnace  with  temperature  controlled  to 
an  accuracy  of  +50,  -25  K.  We  measured  gas  pressure  with  transducers  accurate  to 
better  than  1%  of  the  indication  and  computed  amagat  density  units  using  virial 
coefficients  given  in  Ref.  17.  The  gases  were  mixed  in  the  sample  chambers  while 
monitoring  partial  pressures  and  the  compositions  were  verified  to  within  ±2%  by 
mass  spectrometry. 

The  most  compelling  experimental  evidence  for  a  radiator-speed-dependent 
mechanism  is  summarized  as  follows:  We  observed  non-Lorentzian  line  profiles  only 
for  perturbers  with  masses  much  greater  than  that  of  the  H2  radiator  (for  Ar,  Kr,  and 
N2,  but  not  for  He  or  H2).  For  such  systems,  a  radiator’s  absolute  speed  largely 
determines  its  relative  speed  in  collisions.^  Strong  variations  in  the  collisional  shift 
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with  relative  speed  are  evidenced  by  strong  temperature  dependences  observed  in  the 
density-shift  coefficients.  Thus,  radiators  with  different  speeds  will  be  characterized 
by  different  vibrational  transition  frequencies,  giving  rise  to  an  inhomogeneity.  The 
observed  temperature  dependence  of  the  asymmetry  further  supports  the  speed- 
dependent  mechanism.  Finally,  we  observe  rapid  narrowing  of  the  asymmetry  with 
increasing  H2  partial  pressure,  which  we  believe  indicates  the  important  role  of 
speed-changing  collisions  in  determining  the  line  profile.  Only  results  that  critically 
define  the  line-profile  model  will  be  presented  here. 

Experimental  values  of  the  H2  (2(1)  bne  shifts  for  Ar  and  H2  perturbers  at 
temperatures  from  295  K  to  1000  K  are  listed  in  Table  I.  The  Ar  line-shift  coefficients 
were  determined  by  measuring  Q{1)  spectra  in  1:1  mixtures  of  H2  and  Ar  at  densities 
from  14  to  46  amagat.  We  fit  the  data  to  Lorentzian  profiles  (a  good  approximation  for 
this  mixture  composition)  to  obtain  observed  shifts,  then  corrected  these  for  self-shift 
contributions  using  the  pure-H2  coefficients  measured  by  Rahn  and  Rosasco.*® 

We  modeled  the  2(1)  spectral  profile,  including  the  line  shift,  in  H2-Ar 
mixtures  by  assuming  that  the  collisional  shift  depends  on  the  relative  speed  between 
H2  and  the  perturber.  The  speed-changing  rate  of  H2  resulting  from  H2-Ar  collisions 
was  taken  to  be  comparable  to  the  rate  of  vibrational  dephasing  by  Ar,  but  for  H2-H2 
collisions  it  was  assumed  to  be  much  faster  than  the  dephasing  by  H2.  Also,  we 
assumed  the  following  empirical  form  for  the  2(1)  shift  cross-section  averaged  over 
impact  parameter; 

0(Urel)  =aVrcl  +  ^  (D 

^rcl 

where  is  the  relative  speed  between  radiator  and  perturber,  and  a,  b,  and  c  are 
adjustable  parameters.  Here,  o  has  units  of  cm*  Vamagat  per  cm/s.  The  radiator- 
speed-dependent  shift  coefficient  was  obtained  from: 

^(Urad*  ~  j /m(''ra<f^''rel »  ^ )^rel  ^(^rel )  ^''rel > 
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in  cm'Vamagat,  where  Vrei  is  the  relative  velocity.  Note  that  if /^(v,  T)  is  the 
Boltzmann  distribution  of  perturber  velocities,  /m(vrad+Vrel.  ”0  gives  the  distribution 
of  relative  velocities  seen  by  a  radiator  of  velocity  Vrad-  For  comparison  with  the 
observed  line  shifts,  we  integrated  5(t)rad.  F)  over  the  Boltzmann  distribution  of 
radiator  velocities:  6(T )  =  //jVfCvrad  ,T)S  (t^rad  .  T  )  dvrad- 

Two  sets  of  parameters  a,  b,  and  c  in  Eq.  (1)  were  adjusted'^  to  obtain 
agreement  between  8(J  )  and  the  line-shift  measurements  for  each  collider  pair  (H2- 
Ar  and  H2-H2;  see  Table  I).  The  resulting  £T(t>rel)  for  H2-H2  is  in  excellent 
agreement  (better  than  ±1  from  0.(X)5  to  0.1  eV)  with  the  collision-energy- 
dependent  shift  cross-section  reported  by  Kelley  and  Bragg.^ 

We  now  show  that  the  line  profiles  observed  for  the  H2-Ar  system  can  be 
modeled  using  the  speed-dependent  shift  cross-section  (Eq.  1).  Examples  of  the 
non-Lorentzian  nature  of  these  profiles  are  shown  in  Figs,  la  and  Ic,  where  spectra  of 
Q(l)  from  mixtures  of  2(mole)%  H2  in  Ar  at  total  densities  of  13.9  and  13.7  amagat, 
and  temperatures  of  295  K  and  1000  K,  respectively,  are  plotted.  The  solid  curves  are 
best-fit  Lorentzian  profiles  and  associated  fit  residuals.  Note  that  the  residuals 
increase  substantially  at  lOQQ  K,  but  maintain  a  similar  shape.  Similar  effects  were 
observed  for  Q(0)  and  12(2-4). 

An  inhomogeneous  line  profile  is  computed  using  5(Ura(j,  'O  from  Eq.  (2)  in  the 


expression  for  a  speed-dependent  Lorentzian  profiler^-!® 

r  1  _  if  /A/''rad>F) 


=  :r-S(0),  y,  T)  +  c.c.  s 


I  y(T  )p  -  i[co  +  5(Urad .  T  )p] 


dVrad+c.c.  (3) 


Here,  y  is  the  coefficient  for  the  dephasing  width  (HWHM),  co  is  the  detuning  from  the 


zero-density  Raman  frequency,  and  p  is  the  density  in  amagat.  This  expression  is 
similar  to  that  for  a  speed-dependent  Voigt  profile,^*  except  that  we  have  assumed 


that  )<Urad-  F)  -  y(T)  and  that  the  Doppler-shift  term  (k  Vrad)  is  much  smaller  than 
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yp  due  to  motional  narrowing.^O  The  former  assumption  is  based  only  on  the  fact  that 
no  speed-dependence  of  the  broadening  was  required  to  describe  our  results. 

The  inhomogeneous  profile  given  by  Eq.  (3),  which  neglects  speed-changing 
collisions,  is  generally  in  poor  agreement  with  experimental  asymmetric  spectra  in 
that  it  overestimates  linewidth  and  asymmetry.  However,  if  speed-changing 
collisions  are  present  and  occur  more  frequently  than  dephasing  collisions,  scattering 
amplitudes  associated  with  different  speed-groups  can  coherently  interfere,  resulting 
in  spectral  narrowing.  This  narrowing  and  symmetrization  are  analogous  to  that 
produced  in  Doppler  profiles  by  velocity-changing  collisions,  and  can  significantly  alter 
the  shape  of  the  speed-dependent  profile.  We  included  speed-changing  collisions 
using  a  hard-collision  approximation  analogous  to  that  for  Dicke  narrowing.^i 


resulting  in  a  final  expression  for  the  line  shape: 

,/  X  1  S{03,  y+  V,  r ) 

/(tU)  =  —  - - - -  +  C.C., 

2ff  1  -  vpSiO),  Y+  v,T) 


(4) 


where  v  is  the  frequency  of  speed-changing  collisions.  Note  that  as  v  0, 

1(0))  -4  /inhom(<y).  but  as  v  becomes  larger  than  both  yand  the  width  of  /inhom.  Hci>) 
narrows  to  a  Lorentzian  profile  of  half-width  y  Note  that,  unlike  motional  narrowing 
of  the  Doppler  effect,  varying  the  total  density  alone  does  not  affect  the  degree  of 
narrowing,  since  the  width  of  /inhom  and  vp  both  scale  linearly  with  density. 

To  calculate  the  model  spectra  in  Fig.  lb  and  Id  (2%  H2  in  Ar),  we  used 
Eq.  (4)  and  varied  y,  v,  and  a  frequency  offset  for  best  fit.  (The  latter  is  necessary 
because  the  model  can  vary  from  observed  transition  frequencies  by  up  to  0.5  x  10’3 
cm'Vamagat)  The  resulting  yand  v  were  respectively  2.4  and  8.3  at  KXX)  K,  and  2.4 
and  6.3  at  295  K,  in  units  of  lO'^  cm'Vamagat.  The  temperature  dependence  of  v  (the 
speed-changing  frequency)  is  found  to  be  comparable  to  that  of  P,  the  frequency  of 
velocity-changing  collisions.^O  However,  the  magnitude  of  v  is  only  ~8%  of  p, 
implying  that  the  H2  speed  is  conserved  even  after  several  velocity-changing 
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collisions  with  Ar.  The  inefficiency  of  Ar  in  changing  the  speed  of  H2  is  presumably 
due  to  the  large  mismatch  in  mass  and  average  speed.  The  best-fit  values  for  y  (the 
dephasing  frequency)  were  essentially  the  same;  thus,  the  increased  width  observed 
at  1000  K  results  mostly  from  inhomogeneous  broadening.  Remarkably,  we  find  that 
homogeneous  (dephasing)  broadening  produces  only  ~25%  of  the  total  width  in  this 
case. 

The  model  is  also  found  to  successfully  describe  the  observed  narrowing  and 
symmetrization  of  the  profiles  as  H2  concentration  increases  at  constant  temperature. 
We  measured  spectra  from  a  series  of  gas  mixtures  at  constant  total  density  and 
temperature  but  with  varying  composition.  The  anamolies  were  strongest  in  mixtures 
of  0.5%  H2  in  At  (we  did  not  investigate  greater  dilutions),  but  were  nearly 
indiscernible  in  20%  mixtures  at  295  K.  For  50%  H2  at  295  K,  the  profile  was 
Lorentzian  to  within  our  experimental  uncertainty,  ±1%  of  the  peak  intensity.  More¬ 
over,  the  linewidths  of  the  profiles,  as  determined  from  either  best-fit  Lorentzian  or 
spline  functions,  exhibited  remarkable  deviations  from  a  linear  mixing  rule  (i.e.,  a 
dependence  on  the  mole-fraction-weighted  sum  of  pcrturber  broadening  coefficients). 
This  effect  appears  as  a  curvature  in  linewidths  plotted  against  H2  mole  fraction 
(symbols  in  Fig.  2b).  We  have  observed  similar  nonlinear  dependences  in  coherent 
anti-Stokes  Raman  spectroscopy  (CARS)  linewidths  of  H2  in  Kr,  and  in  IRS 
linewidths  of  D2  in  Ar. 

The  line  profiles  corresponding  to  the  IRS  data  in  Fig.  2b  were  analyzed  by 
varying  the  broadening  and  narrowing  coefficients  (yand  vin  Eqs.  3  and  4)  for  best  fit 
to  each  profile  (described  below),  and  by  employing  a  mole-fxaction-weighted  sum  of 
the  previously  determined  H2  and  Ar  shift  functions  l5(Urad.  T)]  for  each  mixture.  The 
results  for  295  K  arc  shown  in  Fig.  2a.  Using  the  linearly  varying  coefficients  shown 
by  the  solid  lines  in  Fig.  2a,  the  model  accurately  describes  the  dependence  of 
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linewidth  on  H2  mojc  fraction  (Fig.  2b)  and  is  also  in  excellent  agreement  with  the 
observed  spectral  profiles  (not  shown). 

Some  of  the  values  for  yand  vindicated  by  data  points  in  Fig.  2a  could  not  be 
obtained  by  unconstrained  fitting.  At  H2  concentrations  of  50%  and  higher  the 
Lorentzian  nature  of  the  experimental  profiles  provided  poor  sensitivity  to  the  relative 
magnitudes  of  yand  v.  Here,  we  used  a  linear  extrapolation  for  v  which  had  an 
intercept  for  pure  H2  equal  to  the  known  velocity-changing  collision  rate,  (5?-^  (This 
line,  shown  in  Fig.  2a,  was  found  to  be  consistent  with  the  v-values  for  low  H2 
concentrations.)  The  associated  coefficients  for  y  were  obtained  by  fitting  with  fixed  v. 
Thus  using  v  =  for  pure  H2,  the  model  predicts  that  -15%  of  the  collision  width  is 
due  to  inhomogeneous  broadening.  For  H2-Ar  collisions,  about  half  of  the  width  is 
found  to  be  inhomogeneous,  and  v  is  predicted  to  be  -10  times  lower  than  the  pure-H2 
value,  corresponding  to  only  -5%  of  p  for  H2-Ar.20 

Similarly  consistent  results  were  obtained  at  1000  K,  although  only  three 
mixtures  were  investigated  (2%,  50%,  and  100%  H2).  For  pure  H2,  80(±15)%  of  the 
total  linewidth  was  determined  to  be  homogeneous,  and  v  was  equal  to  60(±20)%  of  P 
computed  for  1000  K.20 

We  also  investigated  the  dependence  of  the  line  profiles  on  total  density, 
primarily  in  2%  H2-Ar  mixtures  at  -295  K.  Recall  that  the  model  [Eq.  (4)]  predicts  an 
invariance  of  the  line  shape  with  respect  to  density.  That  is,  if  we  normalize  the 
detuning  by  density  (<u  =  <u'p),  then  p  enters  the  line  profile  only  as  a  proportionality 
constant.  Consistent  with  this  property,  we  found  that  to  within  our  uncertainties  the 
best-fit  values  for  yand  v  were  independent  of  density  over  the  range  2.8-14  amagat. 
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Table  I 


Collisional  line-shift  coefficients  S(X )  for  thcQ(l)  transition  of  H2  at  different 
temperatures,  measured  with  IRS,  and  compared  to  fitted  model  results.  The 
coefficients  are  expressed  in  10'3  cm'Vamagat;  figures  in  parentheses  indicate 
experimental  uncertainties  in  the  last  digit(s). 


H2-Ar  H2-H2 


Temperature  (K) 

Exp.  S(T ) 

Theo.  S(T ) 

Exp.  S(T  )a 

Theo.  <5(7 ) 

295 

-12.0(1) 

-11.4 

-3.20(3) 

-3.2 

450 

-8.2(1) 

-8.3 

0.70(1) 

-0.05 

750 

-2.7(2) 

-3.0 

5.85(8) 

5.4 

1000 

0.86(40) 

1.0 

9.16(14) 

9.5 

^Ref.  18  except  for  the  295-K  value,  which  is  from  Ref.  5.  The  value  for  750  K  is 
interpolated  between  measurements  at  725  K  and  1000  K.  Due  to  its  simplicity,  the 
uncertainty  in  the  model  is  ~0.5  x  10*3  cm'Vamagat;  thus  the  disagreement  at  450  K  is 
not  considered  significant. 
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1.  Experimental  IRS  spectra  (data  points)  of  the  Q{\)  transition  of  H2  broadened 
by  Ar  at  295  K  and  1000  K  The  H2  mole  fraction  is  2%  and  the  total  densities 


are  13.9  and  13.7  amagat,  respectively.  The  data  are  compared  to  best-fit 
Lorentzian  profiles  and  to  the  model  profile  described  in  the  text.  The  model 
includes  inhomogeneous  broadening  arising  from  radiator-speed-dependent 
collisioncil  shifts. 

2.  (a)  Best-fit  collision-broadening  (7)  and  -narrowing  (t/)  coefficients,  indicated 

by  symbols,  obtained  from  IRS  spectra  of  the  Q(l)  transition  of  H2/Ar  mixtures 
at  295  K  and  7.3  amagat  total  density.  (These  coefficients  are  used  by  the 
proposed  line-profile  model.)  Also  plotted  are  best-fit  straight  lines  through 
the  data. 

(b)  Spectral  linewidths  determined  from  fitting  Lorentzian  profiles  to  the  data 
(symbols)  and  to  the  theoretical  profiles  (curve)  calculated  using  the 
corresponding  parameters  in  (a).  Tlie  CARS  measurements  will  be  described 
fully  elsewhere. 
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I 


Raman  shift  (cm"^) 


Narrowing 


Broadening 


•  IRS 

-  CARS 

—  Model 


